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( 57) ABSTRACT 

A structured singular value (p) analysis method of comput- 
ing flutter margins has robust stability of a linear aeroelastic 
model with uncertainty operators (A). Flight data is used to 
update the uncertainty operators to accurately account for 
errors in the computed model and the observed range of 
aircraft dynamics of the aircraft under test caused by time- 
varying aircraft parameters, nonlinearities, and flight 
anomalies, such as test nonrepeatability. This p-based 
approach computes predict flutter margins that are worst- 
case with respect to the modeling uncertainty for use in 
determining when the aircraft is approaching a flutter con- 
dition and defining an expanded safe flight envelope for the 
aircraft that is accepted with more confidence than tradi- 
tional methods that do not update the analysis algorithm 
with flight data by introducing p as a flutter margin param- 
eter that presents several advantages over tracking damping 
trends as a measure of a tendency to instability from 
available flight data. 

1 Claim, 11 Drawing Sheets 
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ON-LINE p METHOD FOR ROBUST 
FLUTTER PREDICTION IN EXPANDING A 

SAFE FLIGHT ENVELOPE FOR AN 
AIRCRAFT MODEL UNDER FLIGHT TEST 

ORIGIN OF INVENTION 

The invention disclosed herein was made by employees of 
the United States Government and may be manufactured and 
used by or for the Government for governmental purposes 
without the payment of any royalties thereon or therefor. 

HELD OF THE INVENTION 

This invention relates to flight flutter testing, which is the 
process of expanding the envelope that determines a range 
of flight conditions within which an aircraft is safe from 
aeroelastic instabilities. This testing must be done for all 
new and reconfigured aircraft. 

BACKGROUND OF THE INVENTION 

Traditional methods of flight flutter testing analyze sys- 
tem parameters, such as damping levels, that vary' with flight 
condition to monitor aircraft stability. (M. W. Kehoe, A 
Historical Overview of Flight Rutter Testing, NASA-TM- 
4720, October 1995.) A real-time method to estimate the 
damping levels was developed based on a recursive 
prediction-error method. (R. Walker and N. Gupta, Real- 
Time FLutter Analysis,” NASA-CR- 170412, March 1984.) 
This method was extended to improve the estimates by 
considering an extended Kalman filter in the formulation. 
(R. Roy and R. Walker, “Real-Time Flutter Identification,” 
NASA-CR-3933, October 1985.) On-line methods using 
both time-domain and frequency domain characteristics of 
turbulence response data have also been formulated to 
estimate dampings. (C. L. Ruhlin et al., “Evaluation of Four 
Subcritical Response for On-Line Prediction of Rutter 
Onset in Wind Tunnel Tests,” Journal of Aircraft, Vol. 20, 
No. 10, October 1983, pp. 835-840.) These methods moni- 
tored stability at test points, but they were of limited 
usefulness for predicting the onset of aeroelastic flutter 
because damping may be highly nonlinear as flight condi- 
tions vary, so damping trends may indicate stability despite 
proximity to an explosive flutter condition. An alternative 
eigenspace method was formulated based on orthogonality 
between eigenvectors, but this method uses a parameter that, 
similar to damping, indicates stability and may vary non- 
iinearly with flight condition. (D. Afolabi et al., “Rutter 
Prediction Using an Eigenvector Orientation Approach,” 
AlAA Journal , Vol. 36, No. 1, January 1998, pp. 69-74.) 

The concept of predicting the onset of flutter by analyzing 
flight data at subcritical airspeeds has been introduced in 
conjunction with a method for formulating a flutter margin 
envelope. (N. H. Zimmerman and J. T. Weissenburger, 
“Prediction of Flutter Onset Speed Based on Right Testing 
at Subcritical Speeds,” Journal of Aircraft , Vol. 1, No. 4, 
July-August 1964, pp. 190-202.) This method considered 
the interaction of two modes in the flutter mechanism to 
formulate a stability parameter that varied quadratically with 
dynamic pressure. This technique has been extended to 
consider several modes interacting as the flutter mechanism 
in order to demonstrate a prediction method for higher-order 
flutter. (S. J. Price and B. H. K. Lee, “Development and 
Analysis of Flight Rutter Prediction Methods,” AIAA 
Dynamics Specialists Conference (Dallas, Tex.) AIAA-92- 
2101, April 1992, pp. 188-200; S. J. Price and B. H. K. Lee, 
“Evaluation and Extension of the Rutter-Margin Method for 
Right Rutter Prediction,” Journal of Aircraft, Vol. 30, No. 


3, May-June 1993, pp. 395^102; and K. E. Kadmka, “Mul- 
timode Instability Prediction Method,” AIAA Structure , 
Structural Dynamics, and Materials Conference (Orlando, 
Fla.), AIAA-85-0737, April 19185, Volume 2, pp. 453-442.) 

5 These flutter margin testing techniques have been used for 
wind tunnel and flight test programs. (R. M. Bennett, 
“Applications of Zimmerman Rutter-Margin Criterion to a 
Wind-Tunnel Model,” NASA-TM-84545, November 1982 
and H. Katz et al., “F-15 Right Rutter Test Program,” 

10 Flutter Testing Techniques , NASA-SP-415, October 1975, 
pp. 413^-31.) However, the method is of limited applica- 
bility for general flight flutter testing because the assump- 
tions of few modes coupling and the requirement to observe 
those modes may be too restrictive. 

15 Stability parameters were also introduced in determining 
flutter margins that consider an autoregressive moving aver- 
age process to describe the aeroelastic dynamics. One 
parameter was based on determinants from a stability cri- 
terion for discrete-time systems that are excited by random 
20 turbulence. (Y. Matsuzaki and \. Ando, Estimation of 
Rutter Boundary from Random Responses Due to Turbu- 
lence at Subcritical Speeds,” Journal of Aircraft, Vol. 18, 
No. 10, October 1981, pp. 862-868 and Y. Matsuzaki and Y. 
Ando, “Divergence Boundary Prediction from Random 
^5 Responses; NAS's Method.” Journal of Aircraft, Vol. 21, 
No. 6, June 1984, pp. 435-136.) A similar parameter was 
developed by extending the determinant method to consider 
short data segments with assumptions of local stationarity. 

( Y. Matsuzaki and Y. Ando, “Flutter and Divergence Bound- 
50 ary Prediction from Nonstationary Random Responses at 
Increasing Speeds,” AIAA Structures, Structural Dynamics, 
and Materials Conference (Orlando. Fla.), AIAA-85-0691, 
April 1985, Vol. 2, pp. 313-320.) Another extension to this 
method derived a similar stability parameter but relaxed the 
35 requirements for station ariness. ( H. Torii and Y. Matsuzaki, 
“Rutter Boundary Prediction Ba^ed on Nonstationary Data 
Measurement,” Journal of Ain raft. Vol. 34, No. 3, May- 
-June 1997, pp. 427-432. ) These techniques of determining 
flutter margins can be applied to complex systems and 
40 require only turbulence for excitation; however, the flutter 
boundary is computed by extrapolating a nonlinear function 
and may be misleading. 

In view of the foregoing, it is clear that in the past the 
actual flight envelope developed for aircraft operation is 
4 5 essentially determined only by flight testing. The edges of 
the envelope are points where either the aircraft cannot fly 
any faster because of engine limitations or, with a 15% 
margin for error, where the damping trends indicate a flutter 
instability may be near. After flight testing, the envelope thus 
50 empirically determined is used for regular operations. It 
would be desirable to use both the aircraft model computa- 
tions and the test flight data in determining flutter margins in 
order to provide a more expanded and robust flutter margin 
envelope. 


55 


STATEMENT OF THE INVENTION 


An object of this invention is to improve flight test 
efficiency along with maintaining a high level of safety in 
60 the process of producing a robust flutter margin envelope for 
an aircraft model. 

A further object is to provide an on-line method of 
producing a robust flutter margin envelope that is essentially 
model based so that it not only has the desired predictive 
65 nature of a traditional method but also utilizes flight test data 
to obtain the desired accuracy of predictive estimates. This 
combines the strengths of both the traditional p-k method 
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and the new method of on-line estimation of the damping 
method. This new method, referred to hereinafter as a unique 
M method of flight testing for flutter margins is for on-line 
flight test prediction based upon both analysis data of the 
aircraft model and flight test data in order that analysis data 
be updated during the testing procedure for continual cor- 
rection of the aircraft model data. 

What is required is a robust stability approach to formu- 
lating a flutter margin envelope that considers a state-space 
model of the aircraft. This method is based on a formal 
mathematical concept of using a structured singular value, p, 
that guarantees a level of modeling errors to which the 
aircraft is robustly stable as described in a technical paper by 
the present inventors R. C. Lind and M. J. Brenner, “Robust 
Rutter Margin Analysis that Incorporates Right Data," 
NASA-TP, March 1998, and presented orally on Sep. 9, 
1997, the presentation of which has been documented by a 
technical memorandum, NASA r TM-97-206220, titled “A 
Presentation on robust Rutter Margin Analysis and a 
Rutterometer, both of which by this reference are incor- 
porated herein. 

A realistic representation of the modeling errors can be 
formulated by describing differences between predicted 
flight and measured flight data. In one application of this 
invention, the method has been successfully used to com- 
pute flutter margins for an F/A-18 Hornet aircraft modified 
into a Systems Research Aircraft (SRA) at NASA Dryden 
Right Research Center. The robust flutter margins were 
determined form aeroelastic flight data to demonstrate the 
potential errors that may exist in the flutter margins com- 
puted by a traditional p-k analysis. (R. Lind and M. Brenner, 
“Robust Flutter Margins of an F/A-18 Aircraft from 
Aeroelastic Flight Data," Journal of guidance. Control and 
Dynamics , Vol 20, No. 3, May-June 1997, pp. 597-604.) 

In accordance with the present invention, on-line estima- 
tion of flutter margins are computed during a flight test in 
contrast to the prior-art flight flutter testing techniques 
referred to above which actually involve a two-step process 
that first requires a computational analysis of the aircraft 
model to estimate off-line flutter margins and then does a 
flight test of the margins. Although the computational part is a 
often not discussed in conjunction with the actual flight 
testing part, it is in fact a necessary part. 

The basic prior-art procedure that is fairly universal with 
industry and military organizations around the world con- 
sists of two parts: 4 

Part I: Pre-flight estimate of flutter margins 

( 1 ) Generate a computer model that is a ‘best-guess' of the 
aircraft model. 

(2) Compute flutter margins for the aircraft model using 

a well known p-k method. 5 

(3) If margins are too small, redesign the aircraft and 
modify the computer model accordingly. Then repeat 
steps 1 and 2. 

(4) If margins are OK, flight test the aircraft. 

Part II: Right test to determine envelope 5 - 

( 1 ) Take aircraft to test point at flight test condition F, a 
dynamic pressure defined by, for example, altitude and 
airspeed. 

(2) Telemeter flight data to ground control room, unless a 
the following step 3 of this Part II can be carried out by 
computers already mounted in the aircraft. 

(3) Estimate dynamics from flight data such as damping. 

(4) Evaluate levels and trends of dynamics. 

(5) If trends are OK, take aircraft to a flight condition 65 
closer to last estimated aeroelastic flutter condition 
(F~F+A f ) and repeat steps 2, 3, and 4. 


4 

(6) If trends at any flight test condition are not OK, then 
declare that condition F as a flutter margin, a point on 
the edge of the aircraft’s safe flight condition envelope. 
Note that the information from the computation step ( 1 ). Part 
5 I, is not used during the flight test. Part II. Conversely, the 
information from Part II is not used during Part I. This is 
essentially why the computational Part I is not mentioned 
when discussing flight testing Part II. Instead, only the 
computation of margins in Part I are used as a first flutter 
10 mar &in estimate. If the flutter margin computation is good, 
it is assumed there are no obvious problems with the aircraft, 
but even if so, a robust flight envelope has still not been 
determined because the computational results can never be 
completely trusted. This is so because the model can only be 
an approximation to the real aircraft and there are often 
important oddities about the real aircraft resulting in errors 
that are not anticipated. The estimated flutter margins that 
are then adopted as the boundaries of the safe operating 
envelope are only determined by flight testing in Part II 
without any method of correcting for such unanticipated 
20 errors. 

Moreover, there is a dramatic risk associated with that 
flight testing procedure because of the unreliability of look- 
ing at damping trends. Damping does not smoothly change 
as an instability condition is approached: rather, the damping 
25 may often undergo sudden changes. Thus, the damping 
trends may seem good at flight condition F, but there is no 
guarantee that the damping will not sharply change as the 
flight condition is changed to F+A^ and aeroelastic flutter of 
the aircraft may occur. This lack of ability to predict an 
30 instability is the reason for the risk in flight flutter testing in 
this tradition procedure. This results in a greater time and 
cost in flight testing because the envelope of safe operating 
conditions must be expanded more cautiously using very 
small increments of AF. 

In contrast, the present invention comprises a method 
used with a singular value, p,, to predict fight margins from 
flight data at successively higher dynamic pressure condi- 
tions of flight at which an onset of a flutter condition will 
occur. Thus, for given dynamic pressures, such as a given 
to altitude and successively greater airspeeds the point at 
which aeroelastic flutter will occur is predicted, and the 
process is repeated for the full range of dynamic pressures 
at which the aircraft is capable of safely operating, such as 
altitudes and airspeeds. Ideally, each flutter margin predicted 
5 at a given altitude will be the same. To expand a full 
operating envelope, the procedure is repeated at all altitudes 
at which the aircraft is capable of operating. In the p method 
of this invention, the flight data is not used to estimate model 
parameters or to identify a transfer function; rather, the flight 
o data is only used to update the uncertainty description for the 
theoretical model of the aircraft being tested. This approach 
avoids several difficulties in trying to estimate a high-order 
model from flight data that has low signal-to-noise ratio but 
still accounts for time-varying dynamics by updating the 
s uncertainty description to describe changing errors between 
the aircraft and the nominal model. 

The stability parameter, p, that is central to this new flutter 
margin method is essentially linear with changes in flight 
condition. Consequently, instabilities can be accurately pre- 
) dieted. In addition, this method can easily consider realis- 
tically high-order models and does not make assumptions 
about the number or type of modes that interact as the flutter 
mechanism. The process used in a flight flutter test is as 
follows: 

Part I: Pre -flight estimate of flutter margins 
( 1 ) Generate a computer model that is ‘best-guess’ of the 

aircraft model. 
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(2) Computer flutter margins for aircraft model using the 
well known p-k method (optional). 

(3) Compute robust flutter margins for aircraft model 
using the singular value p method. 

(4) If the flight envelope determined by computed flight 
margins is too small under either computation method, 
redesign aircraft and repeat steps 1, 2 and 3. 

(5) If the flight envelope thus computed is OK, flight test 
aircraft. 

Note that step 3 is what distinguishes over the prior art. 

Part II: Flight test to determine actual safe aircraft flight 
envelope 

( 1 ) Take aircraft to a safe test point at flight test condition 
F, a dynamic pressure defined by, for example, altitude 
and airspeed. 

(2) Telemeter flight data to ground control room, unless 
the following step 3 of this Part II can be earned out by 
computers already mounted on the aircraft. 

(3) Compute a predicted flutter point F ;> at a higher 
dynamic pressure using a flutterometer algorithm based 
on both flight data and singular value \i: 

(4) If the difference between the present flight test con- 
dition F, and the predicted flutter point F /? is large, then 
expand the next flight text condition to F=F,+A r and 
repeat steps 1, 2 and 3. 

(5) If the difference is small, then declare the last pre- 
dicted flutter margin ¥ p as a point on the edge of the 

envelope. , 

The algorithm is repeated for all altitudes at which the 
aircraft is capable of operating until a completed and robust 
envelope is developed. 

The main advantage to using this \x algorithm over tradi- 
tional methods is that with each iteration at a new flight 
condition F, a new predicted flutter point is determined, 
whereas in the prior-art damping algorithms that is not the 
case. Moreover, the envelope can be expanded more confi- 
dently and in greater increments of AF because there is not 
so much concern about sudden changes in damping and 
stability. The basic flowchart for the flutterometer algorithm 
is as follows: 

1 . Read flight dada D from test point at flight condition F. 

2. Generate computer model P of aircraft at F. 

3. Generate A to relate P and D. 

4. Compute analysis of (P,A) using the p method. 

5. Compute prediction of flutter margin from p, a struc 
tured singular value from step 2. 

A key point is that the flutterometer algorithm works by 
using both the computed model data and the flight data. This 
is important because using only one set of data is insuffi- 
cient The model is only an approximation, so it cannot be 
trusted completely, and the flight data only indicates the 
current state of the aircraft so it cannot be used alone to 
predict the flutter margins. Using both provides robust flutter 
margin prediction which is most important because it 
enables the safe flight condition envelope to be reliably 
established without the aircraft actually reaching or even 
closely approaching a flutter point. Reliable prediction 
throughout the entire process of developing the flight enve- 
lope is due to the fact that the flutterometer algorithm uses 
flight data to continuously update the structured singular 
value p. 

Another key point is that time-varying changes in the 
aircraft dynamics that give rise to modulating errors are 
taken into account. For example, the weight of a fighter 
plane undergoes significant changes as fuel is consumed 
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during the flight tests, but the flutter margins are updated 
correspondingly. This updating is essentially related to the 
manner in which flight and model data are combined in 
generating an operator A that describes changes in the 
5 model. This is so because the structured singular value p is 
a function of the plant or system P (i.e., aircraft) and is 
defined by the following: 

1 

u{P) = 

10 mino^AjiA e A, det (l - PA) = 0 


where: p(P)=0 if no A exists such that det(I-PA)=0, and A=a 
structured uncertainty operator (modeling errors and 
15 perturbations) which is allowed to lie within a norm- 
bounded set such that 


A={A:[A|L^l}, and 

20 is a set of perturbations such that the size of the perturbations 
HAIL is = 1. This means that for robust stability of the plant, 
no perturbation A greater than size ||A|L= 1 can destabilize 
the plant. Thus, the increase of p is the smallest destabilizing 
perturbation, i.e., an exact measure of stability. 

25 An algorithm is provided for updating the aircraft model 
using the uncertainty operator A during flight testing under 
different conditions, and a procedure that uses the structured 
singular value p to determine if the aircraft system P with its 
set of uncertainty matrix operators A is not invalidated. If so, 

30 the set is adjusted until it is not for robust stability of the 
system. This model validation is then used to generate 
reasonable norm bounds for a sufficient set of uncertainty 
operators scaled such that p is always less than unity. 

Two separate algorithms, local and global, are provided 
35 for using flight data sets to update uncertainty operators 
associated with the aircraft plant model P. In the local 
algorithm, flight data at identical flight conditions is used. 
This is done by independently computing uncertainty 
descriptions, i.e., sets of operators A for models at different 
40 flight conditions, resulting in smaller uncertainty operators 
required for subsonic plants and a less conservative worst- 
case flutter margin. In contrast, the global algorithm uses the 
entire set of flight data at all flight conditions to generate a 
single uncertainty description, i.e., set of operators A for all 
45 normal aircraft models. A disadvantage is possibly more 
conservative flutter margins, but an advantage is that the 
uncertainty description used for computing a flutter margin 
(estimate of distance to a flutter point) is truly a worst case. 

Although this statement of the invention speaks of devel- 
oping the flight envelope for new or reconfigured aircraft by 
testing at elevations with progressively increasing airspeed 
(to progressively increase dynamic pressure) in search for 
robust flutter margins (distance to a flight envelope on a 
graph of altitude versus airspeed), it should be understood 
that it is possible to hold true airspeed constant for each 
flutter margin prediction while progressively increasing 
dynamic pressure by changing altitude or Mach number 
which is a function of air density that depends on altitude 
and air temperature. Consequently, while holding a lower 
true airspeed constant, the computed flutter margin may 
predict a point below sea level. That is a valid prediction, but 
the predicted point is outside the flight envelope to be 
promulgated as the universe of safe operating conditions 
with a conventional 15% margin for error at lower altitudes. 

While the \i method of computing robust flutter margins 
has been developed for promulgating an envelope of safe 
operating conditions, it should be noted that in aircraft 
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having greater computer data processing capability the \i 
method can be used with a display as a cockpit flutterometer 
to provide the pilot with virtually real-time information 
about changes in flutter margin in terms of altitude, airspeed 
and dynamic pressure to indicate how far the aircraft can 
drop before reaching a predicted flutter condition as well as 
damping trend to indicate when the flutter condition is near, 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates diagramatically an information flow- 
chart for traditional and methods to compute flutter 
margins. 

FIG. 2 illustrates a block diagram for the small gain 
theorem. 

FIG. 3 illustrates a block diagram with uncertainty for the 
example system. 

FIG. 4 illustrates a block diagram for robust stability 
analysis of the example system using the small gain theo- 
rem. 

FIG. 5a illustrates a linear fractional transformation F W (P, 
A), and FIG. 5b illustrates a linear fractional transformation 
F,(P,A). 

FIG. 6 a illustrates a family of plants P=P(I+AW) with 
input multiplicative uncertainty, FIG. 6b illustrates a family : 
of plants P=P(I+AW)P with output multiplicative uncer- 
tainty and FIG. 6c illustrates a family of plants P=P+AW 
with additive uncertainty. 

FIG. 7 illustrates a linear fractional transformation system 
robust stability analysis using p. 

FIG. 8 illustrates a linear fractional transformation system 
for nominal stability analysis in the p framework with 
parameterization around perturbation in dynamic pressure. 

FIG. 9 illustrates a linear fractional transformation system 3 
for robust stability analysis in the p framework w ith param- 
eterization around perturbation in dynamic pressure and 
structured uncertainty. 

FIG. 10 illustrates a linear fractional transformation sys- 
tem for robust stability analysis in the p framework with ^ 
parameterization around perturbation in dynamic pressure 
and uncertainty in structural stiffness and damping matrices. 

FIG. 11 illustrates a linear fractional transformation sys- 
tem for robust stability analysis in the p framework with 
parameterization around perturbuation in dynamic pressure ^ 
and uncertainty in A Q and matrices of the state-space 
Q(s) model. 

FIG. 12 illustrates a linear fractional transformation sys- 
tem describing Padd approximation to represent unsteady 
aerodynamic force matrix in the p framework with uncer- 
tainty in lag terms. 

FIG. 13 illustrates a family of plants P=P 0 (I+AW) with 
input multiplicative uncertainty. 

FIG. 14 illustrates transfer functions for example system 55 
with multiplicative uncertainty. 

FIG. 15 illustrates a family of plants P=P 0 +WA with 
additive uncertainty. 

FIG. 16 illustrates transfer functions for example system 
with additive uncertainty. 60 

FIG. 17 illustrates linear fractional transformation system 
with nominal models and associates uncertainty operators. 

FIG. 18 illustrates system responses for hardening spring 
example. 

65 

FIG. 19 illustrates system responses for softening spring 
example. 
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FIG. 20 illustrates system responses for hysteresis 
example. 

FIG. 21 illustrates a linear fractional transformation sys- 
tem for robust stability analysis and model validation with 
5 forcing and measurement signals. 

FIG. 22 information flowchart to generate plant and 
uncertainty operators from a system model and flight data 
with the p method. 

10 FIG. 23 is a flow chart for the over-all p method for robust 
flutter prediction and safe flight envelope expansion of an 
aircraft model under test. 

DETAILED DESCRIPTION OF THE 
I5 INVENTION 

1. INTRODUCTION 

Aeroelastic flutter is a potentially destructive instability 
resulting from an interaction between aerodynamic, inertial, 
>0 anc * structural forces. The stability properties of the aeroelas- 
tic dynamics must be investigated to determine a flight 
envelope that is clear of flutter instabilities for new aircraft 
designs or new configurations of current aircraft. Analytical 
predictions of the onset of flutter must be accurate to reduce 
15 dangers and costs associated with experimental estimation. 
Critical flutter conditions are the points closest to the 
flight envelope at which flutter instabilities occur. This 
concept of closeness is formally defined here as a flutter 
pressure that considers the critical dynamic pressure for a 
0 constant Mach value. Obviously, different flutter measures 
such as a flutter speed can be defined because a unique 
equivalent airspeed is associated with each dynamic pres- 
sure for a given Mach number; however, the following 
definition 1.0. 1 for a flutter pressure will be used to 
5 describe the critical flutter flight conditions. 

Definition l.O.l 

A flutter pressure in the smallest value of dynamic pres- 
sure for which an aircraft at a particular Mach number 
experiences a flutter instability. 

The flutter pressure is used to compute a stability margin, 
or flutter margin, that indicates the distance between the 
flutter pressure and a reference point. A common flutter 
margin, T, considers the difference in dynamic pressure 
between the flutter pressure and a point on the edge of the 
flight envelope. Another common flutter margin, n, consid- 
ers the percentage difference between equivalent airspeeds 
at the flutter condition and a point within the flight envelope. 
Definition 1.0.2 

Aflutter margin relates a measure of distance between the 
flight condition associated with the flutter pressure and a 
reference point. 

The traditional p-k method has been extensively used to 
compute flutter margins for a variety of military and com- 
mercial aircraft. (H. J. Hassig, “An Appropriate True Damp- 
ing Solution of the Flutter Equation by Determinant 
Iteration,” A/AA Journal of Aircraft, Vol. 8, No. 11, Novem- 
ber 1971 , pp. 885 - 889 .) This iterative method uses an 
analytical dynamic model coupled with harmonic motion 
solutions for the unsteady aerodynamic forces. The p-k 
method predicts flutter margins entirely from a theoretical 
model that may not accurately describe the true dynamics of 
the airplane. The resulting flutter margins do not account for 
possible variations between the model and the aircraft. 

The community studying aeroelasticity has identified the 
development of improved methods for characterizing flutter 
margins as a vital research area. Flight testing for envelope 
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expansion incurs dramatic time and cost because stability tr 

margins are not computed with a high level of confidence fi 

using traditional methods. The flutter dynamics often exhibit b 

an explosive behavior that results in a sudden change in n 

stability for a small change in flight conditions. Thus, small 5 si 
errors in predicted margins could have grave consequences e 

for aircraft and crews operating near the flutter conditions. fl 

Several approaches exist for characterizing accurate flut- d 

ter margins using flight data generated by the aircraft. These 
data describe the true dynamics and can be used to generate ]0 ^ 
realistic models and compute confident flutter margins. ^ 

Parameter estimation algorithms have been developed to g 

directly identify an aeroelastic model from the flight data. 

The accuracy of the resulting model can deteriorate as the 
complexity and number of degrees of freedom of the system ^ 
increase and signal-to-noise ratios decrease from optimal 
wind-tunnel conditions to realistic flight levels. Modal fil- ^ 
tering has been introduced in association with parameter 
estimation algorithms to simplify analysis by decoupling the r 

system into a set of first-order responses. This type of 1 

filtering does not guarantee robustness and may not perform 
well for systems with many closely-spaced modal natural 
frequencies that cross and shift as flight conditions change. 

Other approaches toward computing confident flutter mar- 
gins evaluate the robustness of a stability margin with 
respect to changes in the model as an indication of the 25 • 
confidence in that margin. A flutter margin robust to pertur- 
bations to the model is a confident margin because model ] 
inaccuracies do not affect that margin. An algorithm has 
been developed to compute the most critical flutter margin 
with respect to first-order perturbations in a model. This 30 
method considers only parametric perturbations and can be 
computationally expensive. A robust control framework has 
been adopted using a feedback structure to relate the struc- 
tural model and the aerodynamic model. This approach uses 
highly conservative robustness conditions with respect to an 35 
uncertainty structure that may not be physically meaningful. 

A similar approach is adopted allowing unmodeled dynam- 
ics and high-order parametric perturbations based on series 
expansion. Statistical approaches are also considered to 
formulate a flutter probability measure. These approaches 40 
will converge to a robustness indicating using Monte Carlo 
simulations, but the computation time can be prohibitive for 
complex systems. The robustness measures for these per- 
turbation and statistical approaches are suspect because no 
global guarantees can be made as to perturbations not 45 
explicitly considered by the minimization algorithms or the 
Monte Carlo simulations. 

An approach to computing flutter margins that guarantees 
a level of robustness and directly accounts for flight data is 
presented herein. (R. Lind and M. Brenner, Robust Flutter 50 
Margins of an F/A-18 Aircraft from Aeroelastic Right 
Data,” AIAA Journal of Guidance . Control and Dynamics , 

Vol. 20, No. 3, May-June 1997, pp. 597-604.) An aeroelas- 
tic model is formulated in a formal robust stability frame- 
work that uses a set of norm-bounded operators. A, to 55 
describe modeling errors and uncertainty. A multivariable 
robust stability measure known as the structured singular 
value, p, computes flutter pressure that are robust to the 
amount of modeling errors as determined by A. (G. J. Bales 
et al., y.- Analysis and Synthesis Toolbox , Musyn Inc. and The 60 
Math Works Inc., Minneapolis, Minn, and Natick, Mass., 
1995.) A robust flutter margin problem is posed by ques- 
tioning what is the largest increase in dynamic pressure for 
which the plant is stable despite possible modeling errors 
described by A. 65 

Right data are easily incorporated into the analysis pro- 
cedure. The modeling errors are determined by comparing 


transfer functions obtained by flight data with transfer 
functions predicted by the analytical model. The norm 
bound on A is chosen based on these observed errors. A 
model validation condition is used to ensure the A is 
sufficient to account for multiple data sets without being 
excessively conservative. With respect to A, a worst-case 
flutter boundary is computed that directly accounts for flight 
data. 

This method illustrated diagrammatically in RG. 1 is 
inherently different from traditional algorithms based on p-k 
methods or parameter identification and robustness 
approaches. The p method uses information from both an 
analytical system model and flight data; traditional 
approaches use only one of these sources, namely the system 
model in the p-k method. Methods that use only an analyti- 
cal model can be inaccurate, and methods that use only the 
flight data can fail if the data are of poor quality. The p 
method uses the flight data to improve the analytical model 
by adding uncertainty operators. Poor quality flight data will 
merely increase the difficulty of obtaining a reasonable 
uncertainty description resulting in a small A. The robust 
margins will be similar to the nominal margins in this case, 
which makes intuitive sense because any information 
obtained from the data should only enhance the plant model 
and improve the accuracy of the flutter margin. 

The concept of computing robustness in flutter margins 
has been recognized for its importance and has recently been 
termed a state-of-the-art research area in aeroelasticity. 
Informal measures of robustness are not necessarily useful 
1 because the informal measures provide no guarantee for the 
system stability. The p method is based on operator theory 
and provides a well-defined concept of robustness that has a 
clear set of guarantees as to the stability properties of the 
system. 


ROBUST STABILITY 

A common definition of a signal is a Lebsegue measurable 
function that maps the space of real numbers R into R . A 
space of such signals is denoted S. 

Definition 2.1.1 

The space of signals that are Lebesgue measurable func- 
tions is S. 

$={/:R— »R'‘} (l) 

Analog measurements x(t) of physical systems are real 
vector functions of the real parameter t describing time and 
thus are valid members of the space of signals, x(t)eS. 
Values of the time parameter, which are often arbitrarily 
numbered as a distance from some reference point, actually 
extend to positive and negative infinity. Stability for physical 
systems must ensure stability for all values of time. A 
time-domain 2-norm is defined as a measure of size (or 
energy) for time-domain singlas x(t)eS that considers all 
time. 

Definition 2.1.2 

The 2-norm measures the energy of a signal x(t)eS. 


One characteristic of a stable system is that only finite- 
energy output signals are generated in response to finite- 
energy input signals. Signals with finite energy are known as 
“square integrable” because the integral of the square of the 
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signal is finite. The Lebsegue space of square integrable 
signals is defined as ^ 2 (— oo©o). This space is also referred to 
as the infinite -horizon Lebesgue 2 -space to denote that the 
norm uses an integral over infinite time. 

Definition 2.1.3 

The space 9L(-oo oo) consists of square integrable time- 
domain signals. 

S£ 2 ( -<«.«> (3) 

Signals associated with physical systems are only known 
for values of time greater than the time at which measure- 
ments are started. Stability analysis and norm computations 
using these signals cannot use properties of the signal before 
the starting time because no information is known. The 
traditional method of characterizing these signals is to 
assume the signal is identically zero for all times before the 
starting time. The time value at which measurements are 
started can be chosen without loss of generality and is 
usually chosen to be t=0. The space ££ 2 (0,oo) is defined as a 
subset of the infinite -horizon Lebesque 2-space to empha- 
size such signals. They are identically zero for all t<0. 

Definition 2 . 1.4 

The space S£ 2 [ 0 ,o°)c: ^£ 2 (— oo, oo) consists of square inte- 
grable time-domain signals that are identically zero for all 
t< 0 . 


12 

notation by rarely distinguishing between time-domain and 
frequency domain signals except where the context does not 
make it clear. The notations for the 2-norm of domain and 
frequency domain are also not distinguished because the 
5 notations are equivalent, as demonstrated by Parse val’s 
identity. 

An important subspace of is the Hardy space, a,. This 
space contains the complex variable tions that are analytic in 
the open right-half of the complex plane and have finite 
10 2 -norm. 

Definition 2.1.8: The Hardy space, W 2 S 2 2 , consists of the 
following functions. 

X 2 -{f(s):f(s)€^ 2 and f{s) is analytic in Re(s)>0} (8) 

15 

System Plant 

A system P is defined as an operator mapping the space of 
input signals S,„ to the space of output signals S owr . This 
definition implies that for any w € S„, and z=Pw, then zeS 

>o oun 

P:S Jf1 — ^ 

Linear, time-invariant systems defined by state-space equa- 
tions are considered. 

1^)1 \A P B P If x{t) I ( i 0 ) 

[ y(r) J [ C P D P ]| u{t) | 


Lz [0. oo ) = |\U): ||.v(/)|| 2 dr <oo, x{t) = 0 for all t < oj 


(4) 


A similar space ££ 2 (-«>,0) is defined for signals x(t) that 
are assumed to begin at t=-oo and are identically zero for all 
times t=0. Thus, the integral to compute the energy for 
elements of this space considers t=-°<> until t> 0 . 

Frequency-domain signals are often considered in stabil- 
ity analysis but do not fall into the set of signals, S. These 
signals /(jto) are complex-valued functions of the imaginary 
unit j=vry, and the real frequency variable co is expressed in 
rad/sec. The set, S ytl) is defined for frequency-domain sig- 
nals. 

Definition 2.1.5 

The space of frequency domain signals is S /0) . 


and P(jo>>=f '(-jo))} <5j 

A frequency-domain 2-norm is formulated to compute a 
measure of energy. 

Definition 2.1.6 

The 2-norm measures the energy of the signal /(jo))€S ya} . 


H/H; = (^P )' 


( 6 ) 


The signal xUteR"' is the state vector, u(t) eR"' is the input 
,0 vector, and y(t) eR" is the output vector. The state update 
matrix is A P eR' M " s ; B^> eR"*"* determines how the input 
affects the states; C P computes the outputs as a linear 

combination of states; and D P eR"-™- [ s the direct 
feedthrough from inputs to outputs. The operator S={A P , 
35 Bp, Cp, D^} denotes the time-domain state system. 

Linear time -invariant state-space system are commonly 
represented by transfer-function operators. These function, 
P(s), are complex-valued matrices of the complex Laplace 
fransform variable, s. Such a transfer- function matrix exists 
40 if and only if the state-space system is linear and time- 
invariant. 

PCs^Dp+CpCsI^— (II) 

45 Stability must be considered over the infinite -horizon time 
lengths so that the operators used map $£ 2 (- 00 , 00 ) into ^ 
(- 00 , 00 ) Properties of the Fourier transform relating ( 00 , 00 ) 
and S £ 2 simply a state-space system, S: $£., (- 00 , 00 )— 
(-oox,oo), is linear and time-invariant if and only if the 
50 associated transfer-function matrix P is such that y=Pu €££•>. 
This condition leads to consideration of the gain for these 
signals. 


55 


IWb = \\Pu\\ 2 
liwlb ~ iiuii 2 


( 12 ) 


A frequency-domain Lebesgue space,$£ 2 , is defined for 
finite-energy signals. 

Definition 2.1.7 

The space S £ 2 consists of frequency -domain signals with 
finite energy. 

^2={/(j<0):/€S / J!/|( 2 <oo} (7) 

The spaces S£ 2 ( - 00 , 00 ) and S £ 2 are isomorphic Hilbert 
spaces under the appropriate inner produce through the 
Fourier transform, which means the spaces have equivalent 
algebraic properties. TN relationship is used to simplify 


This ration of 2-norms will be finite if the system is stable. 
Properties of the 2-norm are used to derive a condition on 
60 the system transfer-function operator, P. This condition is 
referred to as an induced norm because it results from 
consideration of signal norms associated with the operator. 
The S£ 2 induced norm is defined as the H °°-norm. 

65 ^ >e ^ n ^ on 2.2.1: Define the Tf^-norm for transfer-function 
operators. 

|!Pl:^=sup5(p(j(o)) (13) 
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A space of operators with finite „-norm is denoted as S£„. 
Definition 2.2.2: The space consists of system with 
finite df^-norm. 


a_32(P:|'P|i. 


(14) 


14 

The second condition in theorem 2.3.1 is associated with 
the first condition guaranteeing a well-posed and stable 
system. This uniqueness condition can be understood by 
consideration of the solution y for the loop equations shown 
in FIG. 2. 


10 


Transfer functions of linear time-invariant systems are stable 
if and only if z=Pw and w eH, implies zedi 2 . This impli- 
cation results from the Laplace isomorphism between Sf 2 [0. 
oo) and Jf 2 space. These transfer functions are shown to be 
analytic in the open right-half complex plane with finite 
M^-norm. Define the space 5C to contain these operators. 
Definition 2.2.3: The space consists of transfer func- 15 
tions of stable, linear, time-invariant system with finite 

Jf^-norm. 

Jt ; ={P: P >s analytic in Re(s0>0 and ||P|L<“) ■ 

A subspace 9tW,„ is often defined for rational elements. 
Definition 2.2.4: The space, <= consists of ratio- 

nal elements of 

SCK_={P: Pe'K„. and P is rational] ( lfl) 


v=P(u+WMl-PA) ‘Pu 


(171 


20 


Transfer-function operators of linear, time-invariant state- ^ 
space systems are rational functions of the Laplace trans- 
form variable, s. These transfer functions P if and 

only if P is stable such that no poles lie in the closed 
right-half plane. The space 91 which may appear to be 
mathematical abstraction, is thus shown to have a physical 30 
interpretation. 31^ is merely the operator theory repre- 
sentation of stable, rational, transfer functions. 

Small Gain Theorem 

Stability of a linear time-invariant system is determined 
by location of all poles in the left-half plane. Robust stability 35 
in the 94^ and p frameworks is determined by considering 
an interconnection of stable operators. The basis for deter- 
mining stability of these interconnections of operators is the 
“small gain theorem. M 

The small gain theorem states that a closed-loop feedback 
system of stable operators is internally stable if the loop gain 
of those operators is stable and bounded by unity. Several 
formulations of the small gain theorem are derived for 
various signals and systems. Theorem 2.3.1 presents the 
formulation used for this application. 

Theorem 2.3.1 (Small Gain Theorem): Given the feedback 
interconnection structure of FIG. 2 for stable transfer- 
function operators P, A: ££;>— >££2 ^ AeStH^; if the 

TC^-norm of the loop gain is bounded by unity such that 5Q 

||PA|U<1, then: . 

1. the closed-loop system is well-posed and internally 

stable. 

2. a unique y, w e9L is associated with each u €S£„. 

This small gain theorem is overly restrictive in the sense ^ 
of requiring P, Ae9lTC«. A more general small gain theorem 
is formulated for operators not restricted to lie in the 
subspace P, Ae9t3C; theorem 2.3.1 is a special case of this 
general theorem. The extended operator space in the general 
small gain theorem allows consideration of robustness for 60 
systems composed of nonlinear and time-varying operators. 
The requirement of considering stable, rational, transfer- 
function operators is explicitly stated in the theorem to 
emphasize that the nominal aeroelastic system considered in 
this paper is assumed to be stable and the flutter margin is 65 
associated with a destabilizing perturbation to that nominal 
system. 


The inverse term, (I-PA) 1 has a magnitude of infinity if 
the norm of PA is allowed to be unity. Such a condition 
would allow the norm of signal y to be infinite despite a 
norm-bounded u input signal. Restricting ||PA|L<I ensures 
the inverse term exists and a unique finite-norm y is gener- 
ated in response to a finite-norm u. The issue of well- 
posedness requires this condition to hold at s=°° and is 
automatically considered by the 9-f^-norm. 

Robust Stability 

The small gain theorem can be directly used to analyze 
robust stability of a plant model with respect to a set of 
perturbations/ These perturbations are used to descnbe 
uncertainty in the analytical plant model caused by errors 
and unmodeled dynamics. Usually, the exact value of the 
modeling error is not known, but a norm-bounded, real 
scalar, a>0, can be placed on the size of that error. Define the 
A of norm-bounded operators describing these perturbations 
that affect the plant P through a feedback relationship. 


A={A:||A|i„^a} 


(18) 


The small gain theorem allows consideration of the entire 
set of possible modeling uncertainties as described by all 
AeA. The M^-norm of the loop gain cannot be explicitly 
computed for these systems because an infinite number of 
loop gains PA generated by the A exists. The triangle 
inequality of norms can be used to generate a sufficient 
condition for robust stability of P. 


||PA|L^riUIA!U 


(19) 


A condition for robust stability of the closed-loop system 
can be stated. 

Lemma 2.4.1: The plant P is robustly stable to the set of 
uncertainty perturbations. A, that enter the system as in FIG. 
2 with ||A|L=a for all AeA if 


mL<v« 


(20) 


Lemma 2.4.1 shows a sufficient, but not necessary, con- 
dition fo robust stability. The structured singular value, p, is 
introduced in the next chapter as a less conservative measure 
of robust stability that is sufficient and necessary. 

An excellent illustrative example has previously been 
presented to demonstrate the issue of robust stability. This 
example uses classical arguments to compute a robust sta- 
bility condition for a simple system that is seen to be 
identical to the robust stability condition generated using the 
small gain theorem and lemma 2.4.1. A similar example is 
given below for the feedback interconnection in FIG. 3. 

The single-input and single-output elements in the nomi- 
nal system model of FIG. 3 are p, which represents the plant 
dynamics; a, which represents actuator dynamics; and k, 
which represents a feedback controller. Each of the nominal 
system elements are stable transfer functions contained in 
S&TO A modeling error exists on the output of the actuator 
that is represented by a multiplicative uncertainty operator, 
on the output of the element a. 
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The transfer function from w to z can be computed as 
follows. 

2=<-< l+akpr‘akp)w ( 2jj 

Internal stability of the closed-loop feedback system is 
equivalent to stability of the feedback system shown in FIG. 

4 a with the operator g=-< 1+akpr‘akp. 

Because the operators 6 g eSft are stable, the Nyquist 
criterion determines the closed-loop system is stable if and , 
only if the Nyquist plot of Sg does not encircle the -1 point. 
This stability condition is equivalent to the following norm 
condition. 

supjg(jo»6(ja»|<l ( 22> 

This condition is an K^-norm condition on the loop gain, 
g5. Thus, classical Nyquist arguments derive an H^-norm 
condition that is equivalent to the stability condition imme- 
diately formulated by applying the small gain theorem. 

closed-loop stability— I (23) 

The error in the actuator command is unknown and 
possibly time-varying, so the operator S is used to allow 
consideration of a range of errors. Assume the actuator is 
weighted such that the range of errors is described by the set 
of perturbations, ||S)L*<1. Lemma 2.4.1 is used to generate a 
condition that ensures the system is robustly stable to all 
actuators errors described by 8. 

closed-loop stability*— jig I L<1 (24 » 3C 

3. STRUCTURED SINGULAR VALUE 
Linear Fractional Transformations 
The linear fractional transformation (LFIF) is a common 
framework suitable for robust stability analysis using argu- 35 
ments based on the small gain theorem. An LFT is an 
interconnection of operators arranged in a feedback con- 
figuration. These operators can be constant matrices, time- 
domain state-space systems, or frequency-varying transfer 
functions. Consider a linear operator p e c (0| - M * > * <, ' +, -> that is 40 
partitioned into four elements. 


B P , Cp, Dp) matrices of a state-space realization. The 
transfer function can be written as an upper-loop LFT 
involving S and the Laplace transform variable s where i/s 
over s is an integrator. 


■V = Ap.\ + Bpu \Ap Bp 

y = Cpx + Dpu C P Dp 


Pis) = Dp + C P (sZ - A P r 1 B P (28) 


The LFT, F W (P,A), is defined as the interconnection matrix 
such that the upper loop of P is closed with the operator 
AeC /|X °'. 

Definition 3.1.1: Given p € c <0,4 ^ >x</ ' +, 2> and AeC' 1 * 0 ' define 50 
F W (RA) 35 the upper-loop LFT of P closed with A such that 
y=F w (P,A) u as in FIG. 5a. 

FjRAM^+P^Atl-P^r'P,, (26) 

A similar LFT is defined as F,(P,A) to represent the 55 
interconnection matrix of the lower loop of P closed with an 
operator AeC :K °\ 

Definition 3.1.2: Given p € c (0l+02)x(l,+/:> and AeC'-* 02 , define 
FJPA) as the lower-loop LFT of P closed with A such that 
y=F, (P,A) u as in FIG. 5b. 60 

F,(P^>=P 11 +P l2 A(I-P> 2 dr , P 2) P7) 


The LFT is a useful framework for analyzing complex 
systems with many feedback and series interconnections of 
operators. Property 3.1.3 shows the main property of LFTs 
that will used herein. This property allows complex systems 
5 of several interconnected LFTs to be expressed as an equiva- 
lent single LFT. The operators of the new LFT are block- 
structured with blocks composed of the individual operators 
of the LFTs from the original system. 

Property 3.1.3: Feedback and series interconnections of 

0 LFTs can be formulated as a single LFT. 

This issue of stability for LFT systems is associated with 
the concept of a well-posed interconnection. Stability analy- 
sis based on the small gain theorem given in theorem 2.3.1 
can be used to guarantee the LFT is stable and well-posed. 

1 Structured Uncertainty 

The concept of uncertainty is formulated as a set of 
norm-bounded operators. A, associated with a nominal plant, 
P, through an LFT feedback relationship. A family of plants, 

1 T arises through consideration of F U {P,A) for every AeA. 
The true plant model is assumed to lie within this family of 
plants. 

Modeling the uncertainty as a norm-bounded operator can 
lead to overly conservative models. The uncertainty descrip- 
tion can be made more accurate by including frequency 
information. Formulating a model of a physical system that 
is accurate at low frequencies but less accurate for repre- 
senting the system response at high frequencies is often 
possible. A frequency-varying transfer function, W, is gen- 
erally associated with each uncertainty element to describe 
magnitude and phase uncertainty as it varies with frequency. 

Uncertainty can enter a system model in a linear fractional 
manner in several general ways. Two typical types of 
uncertainty are termed “multiplicative" and “additive" 
uncertainty. Multiplicative uncertainty can be either on the 
input or output of a system. Systems with these types of 
uncertainty are easily described in block diagram form. FIG. 
6ci shows the LFT for a plant with input multiplicative 
uncertainty. FIG. 6 b shows the plant with output multipli- 
cative uncertainty. FTG. 6c shows additive uncertainty. 

Uncertainty can also be associated with specific elements 
of the system. These parametric uncertainties are usually 
associated with a system operator in a feedback relationship. 
The number of input and output signals of the system 
operator is increased to account for the additional feedback 
signals associated with the uncertainty operator. This opera- 
tion can be demonstrated by considering P generated by a 
system with an unknown pole. 


+ 1)(j + x) 


An ejtampte ° f an 'erconnection that is common in A norm-bounded, real, scalar, uncertainty parameter 5 can 
stability analysis is t; representation of a time dependent 65 be introduced to account for the possible variation in pole 
state-space system as requency- varying transfer function. value. The set of plants can be written in the LFT framework 
Uehne S as the constant matrix whose entries are the (A,„ using this uncertainty operator and definition 3 1 1 
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-1 

-1 

r= 

f„ip. fyp = 

s + 2.5 

1 

7775 

i 



(s + 1K.T+ 2.5) 

(,v + 1)(.V + 2.5) 


(30) 


, < \.6eR> 
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if ||P|U<1. This robustness condition can be overly conser- 
vative because it does not account for structure in the 
uncertainty operator. The structured singular value, p, is 
defined as an alternative measure of robustness. 

Definition 3.3.1 

Given the complex transfer-function matrix Pe 3ft 'K„ and 
associated norm-bounded set of uncertainty operators A, 
define p. 


10 


(32) 


v(P) = 


15 


20 


A complex system with several uncertainty operators can 
be expressed as an LFT with a single uncertainty operator 
using property 3.1.3. This operator is structured as a block- 
diagonal operator with each block associated with the indi 
vidual uncertainty operators. Two main types of uncertainty 
blocks exist. A full-block uncertainty is a matrix with 
unknown elements. This type of block is used to describe 
unstructured uncertainty in a group of signals. ^ 
Definition 3.2.1: A full-block uncertainty, AeC , has 
unknown elements A, y for every ie(l,n) and j€(l,m). 

A repeated-scalar block introduces more structure into the 
uncertainty description than a full block does. Only the 
diagonal elements of the matrix contain unknown elements, 
the remaining elements are zero. Furthermore, the diagonal 25 
elements are identical. This type of uncertain ty is used to 
relate inpu-output signal pairs with the same uncertainty 

parameter. „ nx/7 

Definition 3.2.2: A repeated-scalar block uncertainty AeC 
has zero-valued elements except an unknown parameter 8 30 
along the diagonal such that A=8l rt . A scalar block is a 
repeated-scalar block of dimension 1 . 

The single structured uncertainty block used for robust 
stability analysis is formally defined in terms of these 
blocks. Let integers m, n, and p define the number of real 35 
scalar complex scalar, and complex full blocks respectively. 

Define integers R, R„, such that the i"’ repeated -sc alar 

block of real, parametric uncertainty is of dimension R ; by 
R . Define similar integers C,, . . . , to denote the 
dimension of the complex repeated -sc alar blocks. The struc- 
tured uncertainty description A is assumed to be norm- 
bounded and belonging to the following set. 


min{tf(A):det(/ - PA l - o} 
AeA 


Define n=0 if no AeA exists such that det(I-PA>=0. 

The structured singular value is an exact measure of 
robust stability for systems with structured uncertainty. The 
value of p determines the allowable size of uncertainty 
matrices for which the plant is robustly stable as demon- 
strated in theorem 3.3.2. 

Theorem 3.3.2 . 

Given the system in FIG. 7, P is robustly stable with 
respect to the A, w'hich is norm-bounded by real scalar a 
such that |lA|USa for all AeA if and only if (t(P)<l/a. 

The model P is usually internally weighted such that the 
range of modeling errors is described by the uncertainty set 
A, which is norm-bounded by 1 . 


|jA|L= 1 for all AeA 


(33) 


40 


A = {A = diagitfh, I c Ai A„): 

sf e R. sf € C. A( e C (> 


(31) 
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Real parametric uncertainty is allowed to enter the prob- 
lem as scalar or repeated-scalar blocks. Complex uncertainty 
enters the problem as scalar, repeated-scalar, or full blocks. 
Complex uncertainty parameters allow uncertainty in mag- 
nitude and phase to be modeled; uncertainty in physical 
characteristics can be more accurately modeled with real 
parameters. The robustness analysis will be less conserva- 
tive by accounting for this structure to accurately describe 
the model uncertainty. 

Structured Singular Value 

FIG. 7 shows the general framework for robust stability 
analysis. The plant operator P(s) e3t Jf is a stable, rational, 
transfer-function matrix representing the aeroelastic dynam- 
ics A norm-bounded AeSt 9U is described such that A(s)eA 
describes the modeling errors in P through a feedback 
relationship. 

The robustness of P with respect to the A can be deter- 
mined using the small gain theorem as presented in lemma 
2.4.1. This condition guarantees stability for any value AeA 
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Theorem 3.3.3 presents the specific condition for robust 
stability that will be used for unity norm-bounded uncer- 
tainty sets. 

Theorem 3.3.3 

Given the system in FIG. 7. P is robustly stable with 
respect to the A with i(A|,_.^ I for all AeA if and only if 

A value of p< 1 implies no perturbation within A exists that 
will destabilize the feedback system. This condition can also 
be interpreted as saying the true plant dynamics are stable, 
assuming these dynamics lie within the range generated by 
the nominal model dynamics coupled with the set of mod- 
eling errors. 

Obviously, p is dependent on the block structure of A. I ne 
robust stability properties computed by p will only be 
accurate if a realistic uncertainly operator is chosen. The 
structured singular value may be arbitrarily greater when 
computed with respect to an unstructured uncertainty opera- 
tor as compared to a highly structured uncertainty operator. 
Definition 3.3.1 demonstrates the p condition of theorem 
3,3.3 is equivalent to the small gain condition of lemma 
2.4.1 when the uncertainty is unstructured. 

Unfortunately, p is a difficult quantity to compute. Closed- 
form solutions exist to exactly compute p for only a small 
number of block structures for A. Upper and lower bounds 
are used to compute p for generalized uncertainty block 
structures. 

4. ROBUST FLUTTER MARGINS 
Nominal Aeroelastic Model 

Consider the generalized equation of motion for the 
structural response of the aircraft. 

M[\$]$‘gR+Ct\+^+qQ(s)r\^0 k34) 

For a system with n modes, define MeR* tn as^the mass 
matrix, CeR' LV " as the damping matrix, and KeR'"" as the 
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stiffness matrix. Define qeR as a scalar representing the 
dynamic pressure, and Q(S) € C ,Uf as the matrix of unsteady 
aerodynamic forces. This equation is valid for a particular 
Mach number, with a different Q(s) describing the unsteady 
aerodynamic forces at a different Mach number. s 

Values of the aerodynamic force matrix at distinct fre- 
quencies can be derived using finite-element models of the 
aircraft and panel methods for unsteady force calculations. 
This research uses a computer program developed for NASA 
known as STARS. This code solves the subsonic aerody- io 
namic equations using the doublet-lattice method. The 
supersonic forces are generated using a different approach 
known as the constant panel method. 

Formulating a linear time-invariant representation of the 
aerodynamic forces to incorporate them into the robust 15 
stability framework is desired. Pad£ approximations can be 
used to compute a rational function approximation to the 
transfer-function matrix. 


£>, = A 0 + sA | + s^A-* + — - — 4 , + — f a (35) 20 

S+Pi S + 02 


This form is often referred to as Roger's form. The 
equation presented here includes only two lag terms, 
although more terms can be included. The poles of the lag 25 
terms, P,, and P 2 , are restricted to be real and positive to 
maintain system stability. The matrix elements of Roger's 
form can be computed using a least-squares algorithm to fit 
the frequency-varying aerodynamic data. 

The aerodynamic lag terms can be replaced in the for- 30 
mutation with a finite-dimensional state-space system rep- 
resented by a transfer-function matrix using Kamel's 
method. 


20 

and s . An approximation to these forms allows them to fall 
within the framework of the method used, including a 
high-frequency pole in the nonproper term, such as replac- 
ing s with 


5+ lOOOO' 


would not affect the low-frequency region of interest while 
ensuring stable and proper functions. With the 
approximation, the forms of Roger and Karpel can be shown 
to be subsets of this method. 

Standard frequency-domain system identification algo- 
rithms can generate a system with an arbitrarily large 
number of states. This state dimension does not greatly 
affect the computational cost of the robust stability analysis. 
Extending the robustness analysis to controller synthesis, 
however, places an emphasis on limiting the state dimen- 
sion. imiting the number of states in the identification 
process is not directly considered here, although standard 
model reduction techniques can be used on the state-space 
system to lower the state dimension. 

Generating a state-space representation of the aeroelastic 
system, including the state-space form of the unsteady 
aerodynamic forces, is straightforward. Consider the force 
vector, y, generated by the state vector, rj. Define xeR" Q as 
the vector of aerodynamic states. 


v = Q{ s)rj » 


A Q b Q 1 1 .v | 

Cq D q I j T] | 


(38) 


Using x, formulate the aeroelastic differential equation 


&s)^ 0 +sA 1 +s 2 A 2 C q Is/-A q )-'B qS 06) 35 

Standard system identification algorithms, including 
curve-fitting or least-squares approaches, can be used to 
compute the elements in the state-space portion of the 
formulation. The A, matrices are assumed to be known from 40 
the low-frequency aerodynamic force data or from experi- 
mental wind-tunnel data. 


M[\$]$ gH+<y\+Kr)+$( C^x+D^r) )=Mti+Cti-H n\+qC^9 > 

A state-space system is formulated using the generalized 
states, qand q, and the unsteady aerodynamic states, x. The 
state -update matrix is determined by the following three 
differential equations. 


A matrix fraction approach is also formulated to represent 
the aerodynamic forces as a linear time-invariant system. 
This generalized form computes rational matrix polynomials 45 
in a fractional form using a least-squares algorithm. Roger’s 
form and Karpel s form can be shown to be subsets of the 
matrix fractional form. 

The approach taken in this application is to fit the aero- 
dynamic force matrices to a single, finite-dimensional, state- 50 
space system. This form is most similar to Karpel’ s form, 
except the additional A,, matrices are not explicitly 
accounted for in the formulation. 

[ Aq Bq i 

QW = \c e D 0 \ = D ^C Q (sI-A e r'B e 


Given the number of generalized states, n, and the number 
of aerodynamic states, n Q , define B 0 €R"o‘", 

and D (? eR' lc ” as the elements of the state-space 
and system approximating Q(s). 

Fitting the aerodynamic data to a finite-dimensional state- 
space system is equivalent to fitting each term in the matrix 
to a real, rational, proper, transfer function. This method 
seems to contradict the methods of Roger and Karpel, which 
form nonproper transfer functions caused by the terms in s 


n 


o / 

0 

7 

(40) 

n 


+ qO Q ) -M ] C 
Bq 0 

A Q 

n 

X 



Nominal Aeroelastic Model in the Structured 
Singular Value Framework 

The generalized equation of motion for the nominal 
aeroelastic system can be expressed in a form suitable for 
using p analysis to compute a flutter margin. The flutter 
margin is dependent on the flight condition parameters that 
55 result in a flutter instability, and ^ is defined to be the 
smallest perturbation among the A set that causes an insta- 
bility. The obvious approach to formulating flutter analysis 
in the ft framework is to introduce a perturbation to a flight 
condition parameter and find the smallest perturbation that 
60 causes an instability. 

Essentially, the two subsystems in the nominal aeroelastic 
model are composed of the structural dynamics, involving 
mass, damping, and stiffness matrices; and the unsteady 
aerodynamics scaled by the dynamic pressure. The gener- 
65 alized equation of motion demonstrates the dynamic pres- 
sure linearly affects the dynamics at a constant Mach con- 
dition. Perturbations to dynamic pressure can thus enter the 
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The dimension of the uncertainty block is the dimension 
of the signal z. The state-space equations for P(s) demon- 
strate this dimension is the number of modes in the system 
n. The number of free variables in the p upper-bound 
optimization, and consequently the computational cost J>f P. 
is a function of the uncertainty dimension. In this way , t e 
number of aerodynamic states. n Q , does not directly affect 
the cost of the flutter estimation. The only cost increase 
caused by these additional states is computing the frequency 
response of the state-space matrix, which is generally muc 
lower than the cost of computing p. 

Demonstrating the aeroelastic system formulated in the p 
framework in FIG. 8 is straightforward. The dynamic pres- 
sure parameterization is equivalent to the nominal state- 
space system given in the previous section. Simply compute 
the closed-loop transfer function with no uncertainty, 65=0, 
and the nominal system is recovered. 

Wind-tunnel and ground vibration testing can expenmen- 
tallv determine aerodynamic stiffness and damping matrices 

that are more accurate than the matrices approximated by a 

The signals z and w are introduced into this formulation 2Q ( inl te-element model. These matrices can be incorpora e 
to associate the perturbation in dynamic pressure to the into a nominal state-space model in the h framework. 
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system through a feedback operator in a linear fractional 
manner that is perfectly suited to p analysis. 

Consider an additive perturbation. S^eR. on the nominal 
dynamic pressure, q () . 

s (41) 

Separate terms in the system dynamics that involve 85. 
0=M[\$}$"gK+CT\M K+qD^m+qC 

1 Crl+AT 1 ( K+qoDeW+qoW ‘ C^v ] 

+5q[M~ 1 DqX\+M~ 1 Qpl 
=n+[M - 1 Cf\+M~ 1 ( K+qoDyW+qoM^ ' l+V 

=rj+[M - 1 Ct)+M~ 1 < K+qoD^+qcM' 1 C qX\+w 


10 


15 


(42) 


lO UIV pv. — y ^ * . . , , 

nominal dynamics in a feedback manner. The signal z can be 
generated as an output of the plant because z is a linear 
combination of states. 

DqT\+M~ ] Cq,x (43) 

The signal w is related to z by the dynamic pressure 
perturbation. 

- (44) 

Vt— Oq- 

The state-space aeroelastic model for nominal stability 
analysis in the \i framework is formulated using the state- 
update matrix. This matrix is determined by the dynamics at 
the nominal dynamic pressure, and the additional input and 
output signals that introduce perturbations to the dynamic 
pressure. That perturbation, 85, is not an explicit parameter 
in the state-space model because 8^ only affects the plant 
through a feedback relationship as deten-mned by the sig- 
nals z and w. Define the transfer function P(s) generated by 
state-space matrices such that z=P(s)w. 

q 1 ( 45 ) 

n 
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FIG. 8 shows the feedback interconnection between the 
perturbation in dynamic pressure and the nominal plant 
model parameterized around that perturbation. This inter- 
connection is an LFT, and the small gain condition of lemma 
2 . 4.1 or the p condition of theorem 3 . 3.3 can be directly 
applied to analyze stability with respect to a variation in the 

flight condition q. . , 

Formulating the nominal aeroelastic dynamics in the p 
framework immediately demonstrates the procedure used in 
computing a flutter margin. Traditional flutter analysis algo 
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This procedure considers variations in dynamic P res ^‘ 
for an aeroelastic model at constant Mach number. The 
unsteady aerodynamics are highly nonlinear with variation 
in Mach number, and attempts to model Mach variations in 
an LFT may produce highly conservative flutter margins. 
The p method presented herein is considering flutter margins 
in terms of dynamic pressure as measured along lines ol 
constant Mach. Flutter is a function of the two variables, 
dynamic pressure and Mach, so computing the dynamic 
pressure causing flutter for a dense set of discrete Mac 
values will generate an accurate portrait of the flutter mar- 

Robust Aeroelastic Model in the Structured Singular Value 

Framework „ . . 

A robust aeroelastic model in the p framework can be 
generated by associating uncertainty operators, A, with the 
nominal model and including the parameterization around a 
perturbation in dynamic pressure. These uncertainty opera- 
tors can resemble any of the forms presented in section 3 . 2 , 
including parametric uncertainty and additive and multipli- 
cative representations of dynamic uncertainty. 

Choosing a reasonable uncertainty description is crucial 
for determining a valid robust flutter margin. This choice can 
arise logically from consideration of weaknesses in the 
modeling process, previous experience with aeroelastic 
analysis, and comparison with observed flig."-,- 
The following section 5 on UNCERTAINTY DESCRIP- 
TIONS IN AEROELASTIC MODELS gives a noncompre- 
hensive examination of several obvious uncertainty descrip- 
tions that may be associated with an aeroelastic model. 

jj, e LFT is a valuable framework for formulating the 
robust aeroelastic model so that the model is suitable for p 
analysis. The various system blocks composed of the nomi- 
nal state-space model with associated uncertainties and any 
additional subsystem blocks and their associated uncertain- 
ties can be expressed as a single model and uncertainty 

operato ^ ^ Qws ^ bloclc structU re used for p analysis of the 


COIHUUllllg d nuiwi »*»“*£,***■ - 

rithms such as the p-k method and the p method as applie The structured singular value 

to FIG. 8 are searching for a value of dynamic pressure tha rnmniirp d w ; t h respect to a single, block-diagonal, struc- 

results in a flutter instability. The nominal flutter margin is P .. . n aetnrK'ifinn tn flvtinmic 


rcsuiis ill a Iiuuv J - ., 

question may be posed, which is answered by these meth- 
ods. 

Question 4 . 2.1 
Nominal Flutter Margin 

What is the largest perturbation to dynamic pressure tor 
which the nominal aeroelastic dynamics are stable? 


65 


tured operator that contains the perturbauon to dynamic 
pressure and the structured uncertainty operator along the 
diagonal. The perturbation to dynamic pressure is explicitly 
shown to distinguish 85 from the modeling uncertainty and 
emphasize 85 as a special operator used to describe a range 
of flight conditions. 
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The flutter margin computed for the uncertain system 
(FIG. 9) is a more accurate margin than one computed with 
traditional methods such as p-k. These traditional methods 
address the nominal flutter problem in question 4 . 2 . 1 ; the 
robust flutter margin must consider the effect of the model- 5 
ing uncertainty. The robust flutter margin actually finds the 
smallest perturbation to dynamic pressure for the entire set 
of plants formulated by the interconnection of the nominal 
dynamics and all elements AeA. Question 4.3.1 poses how to 
compute these margins. 10 

Question 4.3.1 

Robust Flutter Margin 

What is the largest perturbation to dynamic pressure for 
which the nominal aeroelastic dynamics are robustly stable 15 
to the entire range of modeling errors as described by the 
norm-bounded A? 


24 


^ ~~ .fHo ^ : unstable} 


”-q 0 +£ is the nominal 
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This question can be answered by computing u for the 
system in FIG. 9. 

Computing a Flutter Margin with the Structured Singular 
Value 

The nominal flutter problem posed by question 4.2.1 and 
the robust flutter problem posed by question 4 . 3.1 can be 
solved as a p computation. The value of p is a sufficient 25 
direct indication of the flutter margin for the nominal sys- 
tem; additional information regarding the norm bound on the 
uncertainty is required to derive the robust flutter margin. 

The nominal aeroelastic model is formulated for stability 
analysis in the g framework in section 4.2 by introducing a 
perturbation, 85 , to dynamic pressure. A nominal flutter 
margin is computed to answer question 4 . 2. 1 by considering 
the smallest value of this perturbation that destabilizes the 
model. The nominal flutter pressure can be directly calcu- 35 
lated by computing p with respect to the perturbation 
operator 8 ^. 

The exact value of g can be analytically formulated to 
compute a nominal flutter pressure because a closed-form 
solution for p with respect to a single, real, scalar operator 40 
is known. This solution is the spectral radius of the 
frequency-varying transfer-function matrix. 


Then fi(P)= 1/8 such that q nultcr n 
flutter pressure and 

r nom =8 

represents the nominal flutter margin answering question 
Proof 

This result is immediately obvious using definition 3 . 3.1 
for p with respect to a scalar uncertainty parameter 5-. The 
H is the inverse of the smallest destabilizing perturbation, 
and d is computed as the smallest value of 8 r that destabi- 
lizes P; thus p(P)=I/5. □ q 

Lemma 4.4.1 indicates a computational strategy to com- 
pute a nominal flutter margin that does not require a search 
over a set of frequency points. The flutter pressure is found 
by increasing values of 8 ^ until an eigenvalue of the state 
matrix of F w (P, 8 r) has a negative real part indicating F (P ,8 
q) is unstable. Algorithm 4.4.2 demonstrates a bisection 
search implementation that efficiently computes upper and 
lower bounds on the minimum destabilizing 8 ^, perturbation 
to within a desired level of accuracy. q 
Algorithm 4.4.2 (nominal flutter margin ); 

Given plant P at nominal dynamic pressure, q affected 

by perturbation 8^ as in FIG. 4.1: 

Define scalars S as bounds on the smallest 

destabilizing 8^. 

Define scalar e >0 for accuracy. 

while <<W ~ <Wr > e) { 

1 

^9 ~ ~ upper + &iowtr ) 


if F u (P, <5- ) has an unstable pole), then 6 U 


else 






Lkw — At 


/J(P) - max p(P(jto) 


The spectral radius of P(jco) is a discontinuous function of 
frequency, so computational algorithms based on searches 
over a finite set of frequency points may not guarantee the 50 
correct computation of robustness values. A small amount of 
complex uncertainty can be added to the real uncertainty that 
allows a continuous p function to be analyzed but introduces 
unrealistic conservatism. An heuristic robustness indicator 
can be substituted for p that considers stability over a 55 
frequency segment but is not considered here. 

A relatively simple approach can be used to compute p 
with respect to a single real parameter by considering the 
destabilizing value of the parameter. A search over the 
parameter space will result in computation of p and the 
desired flutter margin. Lemma 4.4 . 1 presents the principle of 
this approach. 

Lemma 4.4 . 1 

_ Given the plant, P, derived at nominal dynamic pressure, 
q 0 , with a perturbation to dynamic pressure, 8 -, arranged in 
the feedback relationship of FIG. 9 define 8 . q 


Robust flutter maigins that address question 4 . 3.1 cannot 
be computed using algorithm 4.4.2 because an additional 
search over the set of uncertainty operators A must be 
included. These margins must be computed using the aug- 
mented plant P, which includes feedback signals relating the 
perturbation to dynamic pressure and the uncertainty 
description as shown in FIG. 9. Define the block-structured 
set of operators A^, which considers a particular perturba- 
tion 0 ^ and set of operators A describing model 
uncertainty — 


2*- - {a*- : As- = J A .AeA, ||A|U< l| 


A set of plant models, F W (P, A*), exists for each value of 
fq that defines the nominal plant at dynamic pressure q= 8 0 -h 5 
q and variations to that plant caused by the set of uncertain- 
ties A. The robust flutter margin corresponds to the smallest 
perturbation 85 for which an unstable plant, F„(P, A s ), exists 
with \eA & If every member of the set of plants F^(P, A s ) 
is stable, then q=q 0 + 8 ^ is not a flutter pressure and V 
formulated at q is robustly stable to the uncertainty descrip- 
tion A. v 
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The smallest destabilizing perturbation corresponding 
to a robust flutter margin can be computed by a \i compu- 
tation. The [i framework analyzes robustness with respect to 
a single, structured operator, so the operator set used to 
compute a robust flutter margin must contain the set of 
uncertainty operators A along with a range of dynamic 
pressure perturbations. The p will compute the robustness of 
P with respect to this operator set to find the smallest 
destabilizing perturbation to dynamic pressure and the 
smallest destabilizing uncertainty operator. Define the A that 
contains A 6 _ sets for a norm-bounded set of 5^ operators. 


a = <S : A = 


< 5 ^ 0 

0 A 


L Ae A, HAIL ^ 1. 




( 48 ) 


Imposing the norm bound for 5^ operators as || 8 qlL= 1 ma y 
seem overly restrictive because the units of 6 ^ are die same 
as the units of q in the model. This condition_ implies the A 
considers the range of flight conditions q=q±l lbf/fr for 
plants formulated by dynamic pressure in units ot lbf/fr. 
Such a small range of flight conditions is not useful for 
stability analysis unless q 0 is extremely close to the flutter 
pressure. This limitation is avoided by introducing a weight- 
ing function, W-, to the computation of q. 




(49) 


A W->1 allows a large range of flight conditions to be 
considered despite the unity norm-bound constraint on 
This weighting is incorporated into the stability analysis by 
scaling the feedback signals between the ^ operator and the 
plant P to form the scaled plant, P. 


•■a : 


( 50 ) 


A robust flutter margin is computed by analyzing p(P> 
with respect to the A. The robust flutter pressure is deter- 
mined by iterating over scalings W- until the smallest 
pressure q=^ 0 =W- is found for which the P is robustly stable 
to the set of uncertainties A. Theorem 4.4.3 formally dem- 
onstrates this concept. 

Theorem 4.4.3: Given the plant P derived at nominal 
dynamic pressure q 0 with a perturbation to dynamic pressure 
5- and set of uncertainty operators A norm-bounded by one 
arranged in the feedback relationship of FIG. 9, define the 
plant P with real diagonal matrix W- scaling the feedback 
signals relating Sq and P. 


26 

destabilizing perturbation to dynamic pressure_is_at least 6 
which corresponds to dynamic pressure q=qo + Wq^j — 
q +W-. Thus. F is guaranteed to be robustly stable to the 
uncertainty set A for any perturbation to dynamic pressure 
5 less than W-, to W- is a robust flutter margin. 

—>( sufficient) - 

Assume p(P)>l. Define real, scaler a<l such that u (r)- 
si. which implies, from _theorem 33.2, that P is robustly 
stable to all uncertainties A8A with ||A|L«X<1 . Thus, P is not 
10 guaranteed to be stable over the entire range of modeling 
uncertainty defined by the unity norm-bounded set A. so this 
perturbation is not a valid robust flutter margin and does not 

answer question 4.3.1. __ 

Assume pFcl. Define real, a>l such that u(P)= 5 l, which 
is implies, from theorem 3.3.2, that P is_robustly stable to all 
uncertainties AeA with ||A|L<«- Thus- p is robustl y s “ ble *° 
an uncertainty description larger than that defined by the 
unity nonn-bounded set A, so this condition defines a valid 
flutter margin but is not the least conservative robust flutter 
20 margin and does not answer question 43.1. 

The only difference between models P and P results from 
the scaling W- which scales the feedback signals between P 
and 6 -. No external scaling matrix is allowed to affect the 
feedback signals between P and A because A is defined with 
25 a unity norm bound. Computing p of the plant P with an 
additional scaling on the lower-loop signals would consider 
a scaled set of operators A that does not accurately represent 
the uncertainty description. Therefore, P only scales the Oq 
feedback signals. 

30 Theorem 4.4.3 can be modified to compute a nominal 
flutter margin by changing the identity matrix in the scaling 
used to compute P to a zero matrix. This modification 
eliminates the feedback interconnection between the model 
and the uncertainty description A, so p considers only the 
35 nominal dynamics and computes the smallest destabilizing 
perturbation to dynamic pressure and T roi^'nont' 

Including the uncertainty description A ensures the robust 
flutter margin will be no greater than the nominal flutter 
margin. The robust flutter margin considers the model used 
40 to compute the nominal flutter margin that corresponds to 
the uncertainty operator A=0eA and the models that corre- 
spond to the remaining operators AeA. The conservatism in 
the robust flutter margin makes intuitive sense because the 
nominal flutter margin is the worst-case stability boundary 
45 for a single model and the robust flutter margin is the 
worst-case stability boundary for a family of models. 


r mh ^v n< 


(53) 
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Then is the robust flutter pressure if and 

only if p(P)=l. Also, 




(52) 


represents the least conservative robust flutter margin 
answering question 4.3.1. 

Proof 

^ (necessary ) 

The condition p(F)==l implies that the smallest destabi- 
lizing perturbation to P is described by some AeA with ||AjU 
so no destabilizing AeA exists and the smallest positive 


The proof demonstrating the necessary and sufficient 
condition p(P)=l also makes intuitive sense because ih£ 
equality sign is needed to ensure the flutter margin is valid 
without being overly conservative. If p(P)<l, then no AeA 
causes an instability and the flutter margin is too conserva- 
55 t i ve ^ if p(P)<l, then the system is not robust to all modeling 
errors AeA and the flutter margin is not a valid robust flutter 

margin. , 

Theorem 4.4.3 demonstrates a robust flutter margin can be 
computed by determining a scaling matrix W- for which p( 
60 F)=l. Algorithm 4.4.4 implements an iterative approach to 
compute a scaling matrix for which p(P)e(l±e) for some 
desired level of accuracy e. 

Algorithm 4.4.4 (robust flutter margin): 

Given plant P at nominal dynamic pressure q 0 affected by 
65 unity norm-bounded 8 ^ and A as in FIG. 9. 

Define initial weighting W^. 
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Define scalar e>0 for accuracy. 


I K o 

H. , 


white (p(P\ > ] +*)OR ifi[p) < 1 - £){ 


“ </<) + 


The dynamic pressure q 0 defining the flight condition for 
the nominal plant dynamics must be chosen carefully to 
ensure the robust flutter margin computed with algorithm 
4.4.4 is valid. The nature of u and the upper bound is such 
that all norm-bounded operators centered around the origin 
are assumed to be valid perturbations to the plant dynamics, 
so positive and negative perturbations to dynamic pressure 
are considered by the robust stability analysis. Thus, the 
robust flutter margin computed by p<P)=l could correspond 
to either perturbation 6^=1 or &=- 1. 

The actual dynamic pressure at the flutter instability is 
occurs at a negative dynamic pressure that may be unreal- 
istic for classical flutter analysis. The value of the nominal 
dynamic pressure q 0 can slide along the real axis to a large 
value without loss of generality in the p analysis because this 
parameter linearly affects the nominal dynamics. A simple 
approach to ensure that the robust flutter pressure is a 
positive dynamic pressure is demonstrated in algorithm 
4.4.5, which iterates over increasing values of q 0 until the 
scaling associated with the robust flutter margin satisfies W 
4<Qo- 

Algorithm 4.4.5 (robust flutter margin with q 0 iteration): 
Given parameters as in algorithm 4.4.4: 

Given initial value of nominal dynamic pressure q 0 : 

valid jnar gin = FALSE 


while( vai idjmargin == FALSE) j 


compute plant P at nominal dynamic pressure q 0 
compute an d associated W q from algorithm 4.4.4 
if <W f>q 0 l then q () = MW q 
rise valid_margin - TRUE 


The exclusion of low and negative dynamic pressures for 
stability analysis may not be desirable for all applications 
related to aeroelasticity. An example of such an application 
is the analysis of aeroservoelastic dynamics for a high- 
angle-of-attack aircraft that concerns instabilities at low 
dynamic pressures. The procedure for these types of analysis 
chooses a low value of q 0 and finds the scaling corre- 
sponding to 5^=-l, which computes the low dynamic pres- 
sure instability. Algorithm 4.4.5 can be modified for these 
applications. 

The flutter margin computation must allow for an arbi- 
trary structure of operators A, so an upper bound such as the 


function derived in appendix A must be used. Algorithms 
4.4.4 and 4.4.5 can be adapted by replacing the p calculation 
with a p upper-bound calculation to compute flutter margins. 
A search over frequency points is required when using the 
s upper bound (appendix A), so the accuracy of the robust 
flutter margin requires a dense grid of frequencies associated 
with the natural frequencies of the worst-case dynamics to 
be considered. 

A simple approach can be implemented if the natural 
J0 frequency of the unstable dynamics at the nominal flutter 
pressure can be assumed to be similar to the natural fre- 
quency of the unstable dynamics at the robust flutter pres- 
sure. This assumption can be justified if the uncertainty does 
not change the critical flutter mode between the nominal and 
15 robust pressures, which is often true for systems that 
have relatively small levels of uncertainty and clear sepa- 
ration between critical and subcritical modal frequencies. 
Algorithm 4.4.6 presents this approach, which first com- 
putes the frequency of the nominal flutter mode and then 
20 computes a robust flutter margin from the p upper bound 
' evaluated near that frequency. 

h Algorithm 4.4.6 (robust flutter margin with reduced fre- 
n quency grid): 

, Given the system in FIG. 9: 

e 25 Compute frequency co associated with nominal flutter 
e dynamics using algorithm 4.4.2. 

d Define dense frequency grid ft centered around to. 

Compute robust flutter maring from p upper bound evalu- 
s ated at ft using algorithm 4.4.5. 

j 30 A large and dense frequency grid increases the confidence 
that computed robustness measure is actually an upper 
2 bound for p. Algorithm 4.4.6 must be used with caution 

; because the assumptions behind its use may not be satisfied. 

Computing a robust flutter margin with a dense frequency 
1 35 grid for a particular aircraft at several Mach numbers and 
[ comparing the frequencies of the flutter dynamics to 

■ those of the nominal flutter dynamics is recommended. If 
these frequencies are similar, then algorithm 4.4.6 can be 
considered for further analysis at different Mach numbers. 

40 4.5 Properties of the Structured Singular Value as a Flutter 
Margin 

The flutter computation method described herein uses p as 
the worst-case flutter parameter. The structured singular 
value is a much more informative flutter margin than tradi- 
45 tional parameters such as pole location and modal damping, 
so p presents several advantages as a flutter parameter. 

The conservatism introduced by considering the worst- 
case uncertainty perturbation can be interpreted as a measure 
of sensitivity. Robust p values that are significantly different 
so than the nominal flutter margins indicate the plant is highly 
sensitive to modeling errors and changes in flight condition. 

A small perturbation to the system can drastically alter the 
flutter stability properties. Conversely, similarity between 
the robust and nominal flutter margins indicates the aircraft 
55 is not highly sensitive to small perturbations. 

Robustness analysis determines not only the norm of the 
smallest destabilizing perturbation but also the direction. 
This information relates exact perturbations for which the 
system is particularly sensitive. Thus, p can indicate the 
60 worst-cast flutter mechanism, which may naturally extend to 
indicate active and passive control strategies for flutter 
suppression. 

Damping is only truly informative at the point of insta- 
bility because stable damping at a given condition does not 
65 necessarily indicate an increase in dynamic pressure will be 
a stable flight condition. The structured singular value 
computes the smallest destabilizing perturbation, which 
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indicates the nearest flight conditions that will cause a flutter 
instability. In this respect, p is a stability predictor and 
damping is merely a stability indicator. 

These characteristics of p make the worst-case flutter 
algorithm especially valuable for flight test programs. 
Aeroelastic flight data can be measured at a stable flight 
condition and used to evaluate uncertainty operators. Unlike 
damping estimation, the p method does not require the 
aircraft to approach instability for accurate prediction. The p 
can be computed to update the stability margins with respect 
to the new uncertainty levels. The worst-case stability mar- 
gin then indicates what flight conditions can be safely 
considered. 

Safe and efficient expansion of the flight envelope can be 
performed using an on-line implementation of the worst- 
case stability estimation algorithm. Computing p does not 
introduce an excessive computational burden because each 
F/A-18 flutter margin presented herein was derived in less 
than 2 min using standard off-the-shelf hardware and soft- 
ware packages. The predictive nature of \i and the compu- 
tational efficiency allow a flutterometer tool to be developed 
that tracks the flutter margin during a flight test. 

5. UNCERTAINTY DESCRIPTIONS IN 
AEROELASTIC MODELS 
Parametric Uncertainty in Structural Models 

Parametric uncertainty denotes operators that describe 
errors and unmodeled variations in specific elements and 
coefficients in dynamic system equations. Recall the gener- 
alized aeroelastic equation of motion for state vector qeR - 


30 

five. The actual choice of A* and A ( , is problem-dependent 
and can vary with different aircraft and different finite- 
element modeling procedures. 

The uncertainty operators are described in this section as 
5 A -A^R"*" with real elements because the operators 
describe perturbations to the generalized stiffness and damp- 
ing matrices that are usually real. These operators are often 
additionally constrained to be diagonal operators with n 
independent parameters because the generalized stiffness 
10 and damping are often computed as real diagonal matrices. 
The real and diagonal nature of these uncertainties is not 
required for p analysis, so full-block complex uncertainties 
can be used if they better describe the nature of the modeling 

errors. . 

is Also, the weighting functions W^W r eR MK/ ' are presented 
as constant and real matrices because the functions are 
associated with constant and real stiffness and damping 
matrices. These constraints on the nominal plant through a 
feedback relationship. The feedback operation w is a linear 
20 combination of the states of the plant. The weightings can be 
relaxed if the nature of the uncertainty is best described by 
complex and frequency varying weighting functions. 

Substitute the uncertain K and C into the equation of 
motion, including the state-space representation of the 
25 unsteady aerodynamic forces Q(s) presented in section 4.1. 
Introduce a perturbation, 6^, to dynamic pressure, and sepa- 
rate the nominal dynamics from the unknown terms. 


30 


0=M[\$]$'gh+a\(K^D Q n\+qC^x 


M[\$]$gli+CA+Kr\+iQ'' ! )'\=Q 


(54) 


Robust flutter margins computed with the p method are 
strongly affected by the choice of uncertainty descriptions 
associated with these dynamics, so this uncertainty must be 
a reasonable indicator of potential modeling errors. Para- 
metric uncertainty can be directly associated with the struc- 
tural matrices to indicate specific errors in the finite-element 
model. 

Define an operator, that describes additive uncer- 

tainty of a nominal stiffness matrix Kq. Associate a weight- 
ing matrix, W*€R rtX ” with this uncertainty such that a 
stability analysis should consider a range of stiffness matri- 
ces described by all A* with IIAJU^L 


=t)+[M' 1 C^+AT 1 { K 0 +qoDv)T\+cjoM~ ‘<V1 
+SJ AT 1 D q t\+M~ 1 Q>v]+A a -[ AT 1 W^l+ArlAT 'CHVfil 
=t)+[NT 1 Q1+AT '(Ke+qoDcW+qoM - 1 C^x 

=*1+[AT 1 1 OHAT 1 (K 0 +qoO Q )r\+qcM ~ 1 C Q x]+^+w h +w c 

The signals z- and w q are introduced in section 4.2 to 
40 relate the perturbation in dynamic pressure to the nominal 

r . „ , - • T4. „ ntv»n. 


(57) 


plant through a feedback relationship. The feedback opera- 
tion w^zj is used where z- is a linear combination of the 
states of ffie plant. 


45 


- a =a r x D Q y\+vr { CQK 


(58) 


K=K q +W k A k 


(55) Additional signals are introduced to the plant formulation, 

where z„. and w.. are associated with the uncertainty in the 

K A . : „ » ^ .4 the 


silliness mam*, anu ^ c 

uncertainty in the damping matrix. The outputs of the plant 
z and z c are formulated as linear combinations of the states. 


W licit t'K TT a ***'“ . , , 

Parametric uncertainty can also be associated with struc- 50 stiffness matrix, and z r and w c are associated^ with the 
tural elements in a multiplicative relationship. Define an 
operator, Ac*R n ™ that describes multiplicative uncertainty 
of the nominal damping matrix C G . A weighting W^R"™ is 
again associated with the uncertainty such that the antici- 
pated range of damping matrices for robust stability analysis 
is described by all AC with IIA^IU^ b 


55 


Zk=M~'W k i\ 

-^AT'WcCotl 


(59) 


C=C 0 (I+WcA c ) 


(56) 


The choice of additive uncertainty for A* and multiplica- 60 
five uncertainty for A c does not reflect any generalized 
assumptions regarding the proper way to model errors in 
stiffness and damping; rather, each type is included to 
demonstrate the different mathematical derivations. Addi- 
tive and multiplicative operators are common types of 65 
uncertainty models, so demonstrating how these types ot 
uncertainty are associated with a structural model is instruc- 


The feedback mechanism to describe the modeling uncer- 
tainty uses a relationship between these output signals and 
the w* and w r input signals. 


w*=a*Z* 


(60) 


The state-space plant matrix can be formulated using these 
additional input and output signals. 
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( 61 ) 


FIG. 10 shows how the uncertainty operators and pertur- 
bation to dynamic pressure affect this plant formulation in a 
feedback manner 

The formulation does not directly allow uncertainty in 15 
mass to be described by a feedback operator. The LFT and 
fi frameworks require uncertainty operators to affect the 
nominal dynamics in a linear manner, and this requirement 
precludes introducing mass uncertainty. The inverse of the 
mass matrix scales most terms in the state matrix of P(s), 2( 
including terms involving q. Associating a_mass uncertainty 
operator A M with the mass matrix scaling q would introduce 
terms of A^, which is a nonlinear function of uncertainty 
operators and cannot be directly considered by the p method. 

Parametric Uncertainty in Aerodynamic Models 25 

The unsteady aerodynamic forces Q(s) e C can be 
represented as a state-space model with n^ states. 

Qi s )=D^C q { sl-A q ) x B q (6 2) 

Parametric uncertainty can be associated with the matrix 30 
elements of this state-space representation to describe errors. 
These errors can result from several sources in the modeling 
procedure, including computational fluid dynamic algo- 
rithms that determine the frequency-varying forces and the 
system identification algorithms that represent the compu- 35 
tational values as a state -space system. 

Define an operator A a e R" QX/,Q to describe uncertainty in 
the state matrix of Q(s). This operator directly affects a 
nominal A Qo and describes errors and variations in the poles 
of the state-space representation of the unsteady aerody- 40 
namic forces. Include a weighting function W A e R"*” such 
that the range of state matrices to be considered by robust 
stability analysis is described by all A a ^ with ||A a JL^l. 


which type of uncertainty is most suitable to describe errors 
in A q and is problem-dependent. 

Also, defining the uncertainty operators as real and 
weightings that are real and constant is not a requirement. 
This section presents the problem formulation with these 
definitions because associating these types of uncertainties 
and weightings with the constant real A Q and matrices 
makes intuitive sense. Certainly p can also compute robust 
flutter margins with respect to complex frequency- varying 
weighted uncertainties. 

The nominal aeroelastic model in section 4.2 defined the 
vector x as the states associated with the state-space model 
Q(s). The matrices A Q and b Q directly affect the aeroelastic 
dynamics only through the state derivative equation for x 
and do not appear in the other state derivative equations. The 
aeroelastic plant can be formulated in the u framework by 
substituting the uncertain values of A g and B 0 into this state 
derivative equation without considering changes in the 
derivative equations for-the remaining states r) and t| 

^co +w VV t+B eoff +H VV , > 

^ + * ot ,h+%vv +A vvw< 

Qj\+& Al p. A Q+& B ^BQ 

=A &iL+ B qJ\ +k a. q + w b q 

Several signals are introduced to this equation, where z„ 
and w Aq are associated with the uncertainty in A Q and z* Q 
and w. is associated with the uncertainty in B 0 . The signafs 
z a q ana z « e are output from the plant matrix as Unear 
combination of the states. 


Define also an operator A s e R"^ to describe multipli- 
cative uncertainty in a nominal B^ matrix of Q(s) A 
weighting function W fi( € R"**" is associated with this uncer- 
tainty such that the range of possible B 0 matrices is 50 
described by all with IIA^Jj^l. 


{66) 

The feedback mechanism to describe the modeling uncer- 
tainty uses a relationship between these output signals and 
the and w r input signals. 


The choice of additive and multiplicative operators is 
made for reasons similar to those presented above in the first 
subsection of this section 5. One of each type of uncertainty 
is included to demonstrate the derivation procedures of how 
each uncertainty is associated with the nominal aeroelastic 
dynamics in a feedback relationship. The actual choice of 


A *q Z Aq 

Wb <T Ab q Zs q (67) 

The state-space plant matrix can be formulated using 
these additional input and output signals. 
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PIG 11 shows how the uncertainty operator and pertur- 
bation to dynamic pressure affect this plant formulation in a 
feedback manner. 

Associating a A. uncertainty operator with the B Q matrix 
may not seem immediately useful because considering 
errors in the poles determined by the A t , matrix is often 
intuited. This A B uncertainty can be essential to accurately 
describe the modeling errors because errors m terms com- 
mon to both A G and B Q may exist. Such a situation arises for 
certain modeling representations of aerodynamic lags. Con- 
sider a simplified Roger's form Q(s) matrix that uses two 
Pade approximations to represent lag terms. 



Perform a similar derivation for the state equation of x 2 , 

.V 2 =— (^2* 2~ 2 U Q 

=-(po 0 +W^A {i2 ).\ : -(p2 <) - fU V: A li2 )u O 

=-p 2ff v 2 -p 2o M c >-A 0; (Wj ir Y 2 +Hp 2 M c) ) 

20 

The signals z„ and w P| are introduced as plant output and 
25 input signals to relate the uncertainty A B| in a feedback 
manner. The signals z p , and w p are similarly introduced to 
relate the uncertainty A B , in a feedback manner. The state- 
space matrix can be formulated to describe the nomin 
q o (s) with these additional input and output signals. 


The poles of this system are determined entirely by the A G ,0 
matrix, so uncertainty in the poles can be entirely described 
by a A* operator associated with the A Q matrix. A similar 
uncertainty A B should also be associated with the B G matrix 
in this case because the p„ and p 2 terms appear in both A G 
and Allowing variation in A Q but not in & Q will 
introduce unwanted zeros to the system, so the proper way 
to model pole uncertainty for this formulation is to include 
both A. and A B operators. 

The f’add approximation appears often in aeroelastic 
models, so demonstrating the LFT formulation of Q(s), 
which includes the uncertainties in poles p, and p 2 m a 
feedback relationship, may be useful. The uncertainties in 
CHs) can be developed distinctly from the structural uncer- 
tainties because LFT operations allow the structural and 
aerodynamic models to be combined into a single plant 
model with a structured uncertainty description. 

Define real scalar operators A B] and A B to describe 
uncertainty associated with nominal values of the poles Pi, 
and p, . Real scalar weightings W Pi and W p , normalize the 
uncertainty such that the range of poles to be considered by 50 
robust stability analysis is described by all A B , with 
||A b JLS 1 and all A p , with HApJU= 1- 



Dynamic Uncertainty 

Dynamic uncertainty operators are often associated with 
aeroelastic models to account for modeling errors that are 
not efficiently described by parametric uncertainty. Unmod- 
eled dynamics and inaccurate mode shapes are examples ot 
modeling errors that can be described with less conservatism 
by dynamic uncertainties than with parametric uncertainties. 
These dynamic uncertainties are typically complex in order 
to represent errors in both magnitude and phase of signals 
Consider a system P having two modes with natural 
frequencies at 4 rad/sec and 30 rad/sec. 

/ 16 v 900 \ (74) 

\ s 2 + 0.4$ + 1 6 A 5~ + 0.6 + 900 / 


Pl=Pl 0 +W PlAi| 


Define a nominal model P 0 , vyhich will be used for 
(70) 55 stability analysis of the system but does not model the 
high-frequency mode of R 


Define states x, and x 2 of Q(s), and consider an input 
signal u 0 that generates the output signal y e by the rela- 
tionship y G =Qu G . Substitute the uncertain P, into the state 
equation of 

Ai=-IV,-l V(? 

=-d3i 0 +VV^iA^ilx 1 -(^ l0 +Wp l Ap | >u e 

— P 1 o A I-P 1 o w cr A Pi ( W Pr v > +VV Pi "<?* 

=~fii <r v i-Pi 0 U C»" A t i i C Pi 


P j ±—) 

+ 0.45 + 16; 

60 

The large difference in natural frequency between the 
high-and low-frequency modes of P precludes parametric 
uncertainty associated with the low-frequency mode of P Q 
65 from being a reasonable description of the modeling errors 
in P 0 - The magnitude of any parametric uncertainty associ- 
ated with the low-frequency mode would need to be 
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extremely large to account for the unmodeled high- 
frequency dynamics, so the stability analysis would be 
relatively meaningless because the large uncertainty would 
imply the plant is not accurate at any frequencies. 

A multiplicative uncertainty operator A e C, can be used s 
to describe the high-frequency modeling error without intro- 
ducing the excessive conservatism resulting from a para- 
metric uncertainty description. Associate a complex, scalar 
weighting function W(s) e C with this uncertainty to reflect 
the frequencyvarying levels of modeling errors. 10 


3*={A.+»A||A|Lgl} 


36 


W = 100 


I + 0.5 
s + 500 


(76) 


The set of plants IP used for robust stability analysis is 
formulated to account for the range of dynamics as 
described by the norm-bounded multiplicative uncertainty 


15 


20 


5>={Po(/+tVA):||A|L^l} (7?) 

FIG. 13 shows the block diagram for robust stability 
analysis of IP and FIG. 14 shows the magnitude of the 
transfer function from input to output for P and P 0 , and the 
maximum of [P 0 (I+W)] which is an upper bound for the 25 
output of P at each frequency. The multiplicative uncertainty 
is able to bound the modeling error at the frequency mode 
without introducing excessive conservatism from large 
uncertainty associated with the low frequency mode. 

Dynamic additive operators may also be required in the 30 
uncertainty description to account for errors that are not 
efficiently described by either multiplicative and parametric 
uncertainties. Modeling enrors associated with a zero of the 
system dynamics are an example of an error that is best 
described by additive uncertainty. Multiplicative operators 35 
are not useful in this case because P o (j<o)=0 at the frequency 
to associated with the zero of nominal model and, 
correspondingly, every member of the set of plants P 0 (jw) 
(l+WAM) at this frequency. Additive uncertainty allows the 
system output for some member of the set of plants to be 40 
nonzero even at frequencies of the zeros of the nominal 
plant. 

Consider a plant P with several poles and zeros. 


D _ 16f r + 0,48s + 64 )/ 900 ^ 

64(^ + 0.4&y + 16 Jls 2 + 0.6?+ 900 J 


(78) 
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Assume the nominal plant P 0 of this plant is similar to P 
and has the correct number of poles and zeros, but the 50 
coefficients of the system equations are incorrect. 


(81 ) 

FIG. 15 shows the block diagram for robust stability 
analysis of IP. 

FIG. 16 shows the magnitude of the transfer function 
from input to output for P and P 0 , and the maximum 
magnitude of |P 0 +WA| at each frequency. The additive 
uncertainty bounds the modeling error at each frequency, 
including the frequencies near the zero of the nominal plant! 
because the output of P is bounded above by the maximum 
magnitude of the members of the set IP. 

These multiplicative and additive uncertainties are par- 
ticularly important when comparing an analytical transfer 
function with experimental transfer functions derived from 
flight data. Analytical models are often computed for a low 
range of frequencies because the high frequencies add 
complexity to the model but do not always affect the stability 
margins of the aircraft. The experimental data may indicate 
a high frequency mode that is not included in the analytical 
model, so a frequency-weighted dynamic multiplicative 
uncertainty can be associated with the model. 

The issue of mode shape uncertainty is often encountered 
when comparing low-frequency-predicted dynamical 
responses with flight data because sensor measurements are 
directly affected by the mode shapes. Both multiplicative 
and additive uncertainties may be required to accurately 
model mode shape errors and account for inaccurate 
response levels (which may be higher than predicted at some 
frequencies but lower than predicted at others) and inaccu- 
rate frequencies associated with poles and zeros. 

Uncertainty Associated with Nonlinearities 
The ft framework described uses linear operators to rep- 
resent dynamical models and associated uncertainties but 
does not admit nonlinear operators. The [t framework is 
useful for analyzing aircraft stability despite the constraints 
of linearity because physical systems, which are always 
nonlinear, can often be approximated by linear models to a 
high degree of accuracy. A classic example of this situation 
notes the linearized dynamics are often an acceptable rep- 
resentation of an aircraft operating near trim flight 
conditions, so linear models work well in practice for control 
synthesis and stability analysis. 

Nonlinear dynamics cannot always be accurately 
described by a linearized dynamics model, so the stability 
analysis should consider the effects of these nonlinearities. 
The ft framework can associate linear uncertainty operators 
with linear models to describe the enrors resulting from some 
types of unmodeled nonlinear dynamics. The uncertainty 
does not actually represent the nonlinearity; rather, the 
uncertainty allows the linear system responses to vary with 
sufficient magnitude to bound a range of nonlinear system 
responses. 


r + 0.36s + 36 

900 , 

Os- +0.4s+ 16 

V s 2 + 0,6s + 900 / 
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Define an additive uncertainty operator A € C that is 
normalized by a complex, scalar function W to reflect the 
frequency-varying levels of the modeling errors. 

60 


H' = 0.i 


s + 50 
s + 5 


(80) 


The set of plants IP used for robust stability analysis is 65 
formulated to account for the range of dynamics as 
described by the norm-bounded additive uncertainty A. 


Separating the nonlinear dynamics that cannot be linear- 
ized from the nonlinear dynamics that can be accurately 
represented by linear models is useful. Separate uncertainty 
descriptions can be formulated for each dynamical block, 
and the resulting operators can be combined using the LFT 
framework to formulate a single, linear, plant model with a 
structured uncertainty description. FIG. 17 shows an 
example LFT system representation for a nominal plant P 0 
and associated uncertainty A, and a nominal linear model N 
and associated uncertainty A^ representing a system element 
with nonlinear dynamics. 

The system shown in FIG. 17 is commonly used to 
describe the coupling between the aeroelastic dynamics and 
actuators affecting an aircraft through control surfaces. 
Actuators can display many types of nonlinear behaviors 
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and should be considered in the aeroelastic stability analysis, 
because pilot and autopilot commands that maintain trim 
during flight ensure the control surfaces are continuously 
moving. The errors in linear models resulting from unmod- 
eled actuator dynamics such as nonlinear stiffness param- 
eters or hysteresis functions can sometimes be described by 
a linear uncertainty' operator. 

Consider the response y of a nonlinear system N that 
models a system that has a nonlinear stiffness corresponding 
to a hardening spring that is valid for the bounded input 
signal u e [-10, 10]. Such a system can represent an element 
of an actuator model or a nonlinear structural model. 


N„=2 


( 88 ) 


Associate an additive weighting operator A„ with N 0 such 
that stability analysis considers the set of plants N. 


AWV U + A v : 114x11°°= 1 


(S 4 )) 


^=M/=2ih0.02m i -K).0082« 1 


(82) 


Define a linear nominal model N 0 such that y-N 0 u 
approximates the response of N. 


N 0 =2 


(83) 


Associate an additive weighting operator A„ with N 0 such 
that stability analysis considers the set of plants N. 


N={NoA n 11A^\} 


(84) 


FIG. 18 shows that the maximum and minimum magni- 
tudes of the responses of the set N are able to bound the 
response of the nonlinear system N for the range u e [10, 10]. 
These bounds are overly conservative throughout this oper- 
ating range, but they achieve the desired goal of describing 
errors in the linear system response resulting from the 
unmodeled nonlinearity. 

A similar procedure can be used to describe error caused 
by an unmodeled nonlinear softening spring. Consider the 
response y of a system represented by N that is valid for an 
input signal u e[-10, 10]. 


\=N=2w-K).005« 3 -0.005m 3 
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Define a linear nominal model N 0 such that y-N 0 u 
approximates the response of N. 


No=2 


(861 


Associate an additive weighting operator A^ with N 0 such 
that stability analysis considers the set of plants N. 


AtAV^:|lAJ^ 1 


(87) 


FIG. 19 shows that the maximum and minimum magni- 
tudes of the responses of the set N are able to bound the 
response of the nonlinear system N for the range u e [-10, 
10]. These bounds are also overly conservative throughout 
this operating range, but they achieve the desired goal of 
describing errors in the linear system response resulting 
from the unmodeled nonlinearity. 

Another nonlinearity that commonly affects aeroelastic 
systems is hysteresis. The response of a hysteretic system 
depends on the trend of the input signal such that an 
increasing input signal generates a certain response, but a 
decreasing input signal generates a different response. Such 
hysteresis dynamics are difficult to express as a simple 
mathematical formula, so for illustrative purposes, assume N 
is a nonlinear hysteresis function whose response y depends 
on the trend of the input signal. Define a linear nominal 
model N 0 whose response y=N 0 u approximates the 
response of the hysteretic N. 


10 FIG. 20 shows that the maximum and minimum magni- 
tudes of the responses of the set N are able to bound the 
response of the nonlinear system N with the hysteresis for 
the range u e [-10, 10]. Again, the bounds resulting from a 
linear uncertainty description are overly conservative 
15 throughout this operating range, but they achieve the desired 
goal of describing errors caused by the unmodeled hysteresis 
nonlinearity. 

Explicitly constraining the operating region of the input 
signal u e [-10, 10] can be important to developing reason- 
20 able uncertainties to describe errors resulting from unmod- 
eled nonlinearities. The errors in the linear model can grow 
excessively large when considering a large range of inputs, 
so the uncertainty magnitude would need to also grow 
excessively large. The conservatism resulting from such a 
25 large uncertainty description may be unacceptable and 
require the input range to be constrained to more reasonable 
limits. 

The uncertainty description, even for a constrained oper- 
ating region, will usually be overly conservative when 
30 describing modeling errors for certain parts of the operating 
region. The uncertainty is able to bound the errors in FIGS. 
19 and 20, but the maximum and minimum bound are 
clearly not optimal. Some amount of conservatism is 
expected when describing errors resulting from unmodeled 
35 nonlinearities because a linear model, whether a single plant 
or a set of plants, will usually not be an accurate represen- 
tation of a nonlinear system that cannot be linearized. 

Also, several types of unmodeled nonlineanties exist that 
are frequently encountered in aircraft systems but are not 
40 easily described by linear uncertainty operators associated 
with the linear models. Examples of these types of nonlin- 
earities include free play, dead-band responses, friction, and 
rate limiting of actuators. 

The robust stability margins computed from g with 
45 respect to the linear uncertainty operators describing unmod- 
eled nonlinear dynamics will always be somewhat suspect, 
because this approach cannot consider stability properties 
unique to nonlinear systems such as bifurcation points and 
limit-cycle behaviors. This approach limits the usefulness of 
50 the p method to systems for which the nonlinearities have 
small effects on the response and do not introduce nonlinear 
instabilities to the critical flutter mechanism. 

Uncertainty Associated With Flight Data 
A theoretical model with an associated uncertainty 
55 description can be an accurate representation of the aeroelas- 
tic dynamics of an aircraft, but responses from that model 
may not identically match flight data. Additional uncertain- 
ties can be associated with the model to describe errors that 
are observed between the predicted responses and the mea- 
60 sured responses from a commanded excitation to the aircraft. 
These uncertainties do not necessarily indicate errors in the 
model; rather, these uncertainties indicate errors in the 
process used to generate aeroelastic responses and measure 
flight data. 

65 One source of error is an incorrect assumption of excita- 
tion force used to generate the predicted and measured 
responses. The measured excitation force associated with 
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the flight data may not correctly account for poor hardware 
performance and nonuniform spectral distribution of the 
force. Also, inexactly phased excitation between multiple 
force mechanisms can excite modes that are not anticipated 
b) a theoretical analysis. A frequency- varying dynamic 
uncertainty can be associated with the force input of the 
analytical model to describe errors in the excitation. 

The phenomenon of nonrepeatibility can cause discrep- 
ancies between predicted and measured responses from 
multiple occurrences of excitation signals. Nonrepeatibility 
affects flight data by introducing slight variations in 
responses, even for data recorded at identical flight condi- 
tions with identical excitation signals. This unexplained 
behavior may result from some unmodeled nonlinear 
dynamic or inexact excitation that is not correctly measured. 
External disturbances such as wind gusts or turbulence can 
introduce an unmodeled dynamic that inconsistently affects 
the aircraft responses. A frequency varying dynamic uncer- 
tainty can be associated with the model to describe nonre- 
peatible data variations. 

Another source of error between predicted an measured 
responses is an incorrect assumption of flight condition. 
Flight data sets are sometimes generated at test points that 
attempt to maintain a constant flight condition to match the 
data sets predicted from a model describing the aeroelastic 
dynamics at that same flight condition. Slight variations in 
flight conditions while the experimental response is mea- 
sured may cause some discrepancy between the predicted 
response and the flight data. Parametric uncertainty associ- 
ated with the unsteady aerodynamic model can be used to 
account for these errors because flight condition variations 
only affect the aerodynamic model and not the structural 
model. 

The model may accurately represent the mode shapes of 
the aircraft but have a poor representation of the sensor 
locations. The responses measured by sensors are inherently 
dependent on sensor location, with respect to mode shapes, 
to determine the magnitude and phase of the signal. Addi- 
tional errors in magnitude and phase are introduced when 
considering transfer-function estimates generated by signals 
that violate assumptions used in computational algorithms. 

A frequency-varying dynamic uncertainty can be associated 
with the output of the plant model to describe errors in both 
mode shape and sensor location predictions. 

The choice of signal -processing algorithms can also intro- 
duce errors between predicted an measured responses. The 
Fourier transform, which is the traditional tool for signal 
processing, assumes several characteristics of the data that 
may be violated with transient-response aeroelastic data. 
Filtering and wavelet-based algorithms may be used to 
reduce errors introduced by signal processing and reduce 
conservatism in the resulting stability margins. Parameter 
uncertainty associated with the modal parameters of the 
linear model may be used to describe some errors in the 
natural frequencies and dampings observed using flight data 
analyzed by incorrect algorithms. Dynamic uncertainty may 
also be required to describe errors introduced by leakage and 
aliasing effects. 
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beyond the true modeling errors and can be overly optimistic 
if the uncertainty does not sufficiently account for the true 
modeling errors. 

Model validation algorithms can be used to indicate if an 
uncertainty description is reasonable with respect to flight 
data. These algorithms consider whether or not a set of data 
measurements could have been generated by a proposed 
model that includes the nominal dynamics and associated 
io uncertainty operators and noise and disturbance signals. 

Uncertainty operators associated with the nominal plant 
model specifically in the LFT framework can be considered 
by validation algorithms. Frequency-domain algorithms are 
generated that consider the model validation problem in the 
15 context °f control design. Tune-domain approaches are also 
developed that can be solved with convex optimization 
algorithms for certain uncertainty structures. 

A model validation procedure has been formulated that 
uses a p condition to determine if an uncertainty model is 
20 invalidated. This procedure uses frequency-domain transfer- 
function data to determine if some perturbation AeA the 
nominal plant could produce the measurements. Consider 
the block diagram for robust stability analysis of systems 
with measurement y and forcing u signals shown in FIG. 21. 
25 The model validation question, as applied to uncertainty 
models, is given in question 6.1.1. 

Question 6.1.1 (uncertainty validation): Is there some 
frequency-varying AeA for FIG. 21 such that F„(P,A) could 
generate the set of observed data y and u? 

30 Define P(s) ec" 0+w>m,+ " v as a stable transfer- function sys- 
tem matrix such that P e & , Tt. Partition this matrix into four 
elements such that P n (s) ec"’™\ P 22 (s) €c"* Wi , and P 21 (s) are 
of appropriate size. P 22 (s) is the nominal-plant transfer 
^ function in the absence of uncertainty. 


p = 


A 1 P\2 
P 2 | P 2 2 


( 90 ) 


Robust stability analysis of this system is determined 
using the transfer function P , , which contains the feedback 
relationship between the plant dynamics and the uncertainty 
operator. The robust stability condition given by theorem 
3.3.3 requires p(P M )<l. 
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The model validation condition uses these elements of P 
and the finite-energy frequency-domain signals y, u eTU 
where the measurements are scalar functions y, u e C. 
Formulate the following two matrices. 
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F| 2 =P] 2 U 

P,2=P. 2 u-y 


MODEL VALIDATION OF UNCERTAINTY 
Model Validation Using the Structured Singular Value 
Robust stability analysis considers the stability of a sys- 
tem subject to a set of modeling errors and perturbations 
represented by a norm-bounded A. A logical question that 
arises is how to reasonably determine this set. This issue is 
important for computing robust flutter margins because p 
can be overly conservative if the uncertainty is excessive 


( 91 ) 

55 The following theorem formulates the model validation test 
using p. 

Theorem 6.1.2: Given measurements y generated by 
inputs p, then the system P with associated uncertainty set, 
A is not invalidated if 
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( 92 ) 

This condition may seem counterintuitive in that the 
desired condition for validation is p greater than one, while 
65 the robust stability condition seeks a value less than one. 
This condition can be explained by considering the follow- 
ing relationship of y=F u (P, A) u as shown in FIG. 6. 1. 
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-continued 


0=[/* 22 «^yl+/ > 2 iA(/-P,,A> '[Pi2 U ^22 +P 2A (! /> n^ ) 


( 93 ) 


Define the plant P={P„, P 12 , P 21 . ^ A value of (t(P)<l 
implies this system is robustly stable to all perturbations 
AeA. This robust stability of \i (P)<1 also implies lhatthe 
loop gain F„(P. A) is not singular for any value AeA. Thus 
u (P)<l would contradict the relationship shown in the 
above equations that requires F„ (P, A) to be singular to 
satisfy the input-to-output relationship of the data. In this 
respect, the model validation test is actually an inverted 
robust stability test. 

The phrase “model validation’' may be misleading. No 
analytical model can ever be truly validated by considering 
a finite set of experimental data. A model may not be 
invalidated by the data, but no guarantee exists that a 
different data set could not be generated from the physical 
system that invalidates the model. The finite sets of mea- 
surement data cannot record the response of the system to an 
infinite number of input signals subject to an infinite number 
of initial conditions. Theorem 6.1.2 reflects this fact by 
explicity stating the condition only determines if the model 
is invalidated by the data. The system is assumed to be 
nearly linear, and the data is assumed to sufficiently repre- 
sent the behavior of the dynamics, so theorem 6.1.2 repre- 
sents a model validation test. 

Validating Norm Bounds for Uncertainty 
The value of p computed using theorem 6.1.2 can be 
interpreted as a measure of how reasonable the uncertainty 
description is. A value of p=2 implies the uncertainty can be 
scaled by 2 before the model is invalidated. Alternatively, a 
value of a=0.5 implies twice as much uncertainty is required 
for the model to not be invalidated. This interpretation of fl, 
as relating to the size of allowable perturbations, emphasizes 
the relationship of the model validation condition to a robust 
stability condition. 

This model validation is used in practice to generate 
reasonable norm bounds for an uncertainty description. The 
following algorithm can be used to determine a sufficient 
level of uncertainty required such that the model is not 
invalidated by multiple data sets. A small scalar col is 
chosen to scale the uncertainty set and increase the amount 
of allowable errors if the size of A is not sufficient. 
Algorithm 6.2.1 (Model Validation): 

Given frequency -domain data sets {y |» y 2 - • • • » y ^ 
{U,, U05 • • • ♦ 

Given frequency -domain transfer-function P with ele- 
ments P u , Pj 2 * ^ 22 : 

Given uncertainty set A with initial norm bound: 

Given update scalar col: 


/>, | = aP 1 1 (equivalent to A - a&) 


valid = FALSE 
while( valid = F ALSE){ 
valid = TRUE 
for / = 1 : n | 

P12 = P\2U i 

? 22 = P 22 M; - yi 

if tt(P u - 'I 

valid = FALSE 

) 


A is the required norm- bounded uncertainty 

The norm bound on the uncertainty set A is not actually 
10 increased with this method as denoted by the parentheses 
around the statement A=aA. Scaling the uncertainty and 
retaining these scalings throughout the procedure would be 
difficult The algorithm is simplified by always considering 
1S the uncertainty is scaled such that u<) is always the desired 
result. This scaling of the uncertainty is actually accom- 
plished by the statement P , , =aP , , which scales the feedback 
signals between the plant and uncertainty operators. The 
norm bound for the uncertainty that is needed to ensure the 
20 uncertainty levels are not invalidated by the data is scaled 
into the plant by the parameter a. Robust flutter margins can 
be computed directly from the new scaled plant using 
algorithm 4.4.6 with a desired pel condition. 

Algorithm 6.2.1 is a straightforward implementation ol 
25 the mo del validation test. The p values are only compared 
with 1 to determine if the model is invalidated. More 
sophisticated algorithms could use the value of p to deter- 
mine the factor a to scale the uncertainty. Also, the values 
of p across frequency could be exploited. Frequencies with 
30 low p value indicate areas where the uncertainty set is least 
conservative; frequencies with high p are overly conserva- 
tive. The scaling a could vary with frequency to reflect this 

information. - 

The situation may arise when the initial value chosen for 
5 the norm bound of the uncertainty set may be overly 
conservative. In this situation, the model validation condi- 
tion will pass for each data set during the first processing of 
the outer loop. The initial norm bound can simply be 
40 decreased by some level consistent with the lowest p value 
computed during the validation checks, and algorithm 
can be run again to compute a less conservative uncertainty 
description that does not invalidate the model. 

Theorem 6.1.2 is only valid for scalar data signals gen- 
45 erated by systems with a single input and single output. 
Multiple data signals can be considered by applying theorem 
6 1 9 to each combination of singleinput and -output sig- 
nals. Algorithm 6.2.1 can still be used by simply including 
outer loops to cycle over the number of input and measure- 

50 ment signals. . 

The model validation algorithm using theorem 6.1— is a 
departure from traditional methods of analyzing flight data 
to assess accuracy of an analytical model. The most widely 
used algorithms for analyzing flight data estimate natural 
55 frequencies and modal damping. The model validation pro- 
cedure using p considers the response dynamics at each 
frequency without explicitly comparing modal properties. 
This procedure enhances the ability of the p method to 
60 analyze flutter stability without requiring damping esti- 
mates. . 

Another interpretation of algorithm 6.2.1 uses the uncer- 
taintv A to bound the magnitude and phase of the possible 
transfer functions generated by the family of plants F„(R )• 
65 The structured singular value ensures the experimental da 
transfer function lies within these analytical bounds at every 
frequency. 
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PROCEDURE FOR THE |i METHOD 

Model Updating 

Generating a model by analyzing flight data is essential 
for computing a confident stability analysis. A nominal 
model generated purely from analytical equations of the 
predicted aircraft dynamics may not accurately describe the 
true aircraft. A model must be generated that accounts for the 
flight data to ensure the predicted dynamics represent the 
true dynamics. 

The most direct method of generating a model from the 
flight data is to identify a system model entirely from the 
data measurements. Many system identification algorithms 
exist that have become standard tools for systems and 
control engineers. Direct application of these methods to 1 
aeroelastic systems rarely produces an accurate model that 
accounts for the dynamics of the aircraft. Aeroelastic 
response data is typically of poor quality relative to ground 
vibration test data because of the low signal-to-noise ratio 
and unobserved dynamics in the response measurements that 2 
may drastically lower the effectiveness of system identifi- 
cation algorithms. 

An alternative method is to use the nominal aircraft 
dynamical model as an initial estimate to model the true 
aircraft. The flight data are then used to update the elements 2 
of this model. Several methods have been devised to update 
an analytical structural model using experimental data. 
Model updating can be performed on the full stress model or 
a subset computed with Guy an reduction. Generally, con- 
sidering the full model is preferable because the reduction * 
may distribute local errors throughout the entire model if an 
orthogonality condition is violated. 

There are two basic methods for updating the full struc- 
tural model using comparisons between experimental and 
predicted data. One method updates the mass and stiffness 35 
matrices of the finite element model. This method suffers 
from lack of physical interpretation of the matrix updates 
and possible numerical conditioning. Another method 
updates specific parameters in the model. This method is 
accurate for small systems but may require an excessive 
computational cost for large systems. 

Aeroelastic models have the additional freedom of updat- 
ing the aerodynamic and the structural elements. Another 
method of updating the linear model in a modem control 
framework has been developed. (K. D. Gondoly, Application 
of Advanced Robustness Analysis to Experimental Flutter , 
Masters of Science in Aeronautics and Astronautics Thesis! 
Massachusetts Institute of Technology, June 1995.) This 
method may be overly conservative for describing 
nonlinearities, and the corresponding stability margins are 
only accurate for flight conditions near the instability. A 
parametric identification algorithm has been developed that 
uses flight data to update specific terms in the aerodynamic 
model through a nonlinear optimization. This method suffers 
with flight data because of unobserved dynamics and the low 
signal-to-noise ratio in the measurements. 

The approach taken in this invention is to update only the 
uncertainty operators of the robust aeroelastic model, using 
the flight data, and leave the nominal dynamics model ( 
unchanged. The model validation condition of theorem 6. 1 .2 
is used with the nominal plant F to generate a reasonable 
uncertainty description A to associate with P. FIG. 22 shows 
the flow of information through the p method. 

Algorithm 6.2. 1 presents a procedure to implement the p 65 
method in a manner corresponding to FIG. 22. The model 
update procedure based on algorithm 6.2.1 actually scales 
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the plant but has the equivalent effect of scaling the norm 
bound of the uncertainty description. The algorithm loops 
over the model validation procedure until an uncertainty 
description is determined that is not invalidated by the flight 
5 data. The final step is to compute a robust flutter margin 
from the scaled plant by using a pel condition as in 
algorithm 4.4.6. 

Algorithm 7.1.1 (robust flutter margins with model 
updating): 

10 Given nominal plant P: 

Given uncertainty set A associated with P: 

Given input excitation data u: 

Given output response measurement data y: 

15 Define W^I as the scaling update for A 

while (F U (P, A) is invalidated by w, y using algorithm 62. 1){ 

A = IV A 


$ "flutter ° b flutter margin computed from algorithm 4 4 6 

flutter 

5 Several advantages exist to using this method as com- 
pared to traditional model updating methods. The typically 
poor quality of flight data, in association with aircraft 
dynamics consisting of many modes, makes updating a 
nominal model difficult. Traditional methods of norm-based 
1 u P dat e algorithms often generate a nonunique set of model 
updates that have no way to determine which update has the 
most logical physical interpretation. The method of updating 
the uncertainty operators based on a worst-case magnitude 
avoids this problem. 

Also, this method can work with flight data of varying 
quality. The updated uncertainty description i highly accu- 
rate if the data show a high signal-to-noise ratio and much 
of the dynamics are observed by the sensois. if the data do 
not have these desired characteristics, however, the method 
can still compute a flutter margin. An uncertainty description 
may be difficult to compute if the data do not indicate the 
aircraft dynamics well, so the model validation procedure 
will not require a large magnitude for the uncertainty opera- 
tors. The robust model in this situation will closely resemble 
the nominal model. In this way, the p method will always 
generate a more accurate flutter margin, and at worst, the 
robust p flutter margin will be equivalent to the nominal 
flutter margin. 

Approaches to Use Right Data 

A flight test generally consists of maneuvers at several test 
points that may be at identical or different flight conditions. 
The entire flight test program will use many flight tests to 
measure response data at test points throughout the flight 
envelope. The model-updating method that generates uncer- 
tainty operators can use the entire set of flight data from the 
different test points. 

Several approaches are formulated to use multiple flight 
data sets to update the uncertainty description associated 
with a nominal plant model. The uncertainty description 
may be different for each approach, and the resulting flutter 
margin will be different for each approach because of the 
dependence of p on the uncertainty set. Two approaches 
discussed here are denoted as local and global. 

A local approach uses flight data from test points at 
identical flight conditions. These data are used to generate an 
uncertainty description for the nominal model at the par- 
ticular flight condition associated with the data. The mag- 
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nitude of the uncertainty operators is chosen such that a 
robust model at the single flight condition is not invalidated 
by any of the flight data sets measured at the same flight 
condition with no consideration of data from other flight 
conditions. 

The local approach shows the benefit of independently 
computing uncertainty descriptions for models at different 
flight conditions. This approach allows each uncertainty 
description to be more accurate because, for example, the 
flight data may indicate much smaller uncertainty operators 
are required for subsonic plant models even though large 
uncertainty operators are required for transonic plant mod- 
els. The resulting worst-case flutter margins will be less 
conservative because less uncertainty is required for the 
model. 

The following algorithm 7.2. 1 outlines the local approach 
to use flight data and compute flutter margins. 

Algorithm 7.2.1 (local approach): 

Define scalar n^as the number of flight conditions. 

Define scalar n, as the number of data sets at flight 
condition i. 

Define vectors {u/, u 2 u" k } as input excitations at 
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The following algorithm 7.2.2 outlines the global 
approach to use flight data and compute flutter margins. 

Algorithm 7.2.2 (global approach): 

Define scalar n^ as the number of flight conditions. 

5 Define scalar n, as the number of data sets at flight 
condition i. 

Define vectors {u,\ u/\ . . . , u/ 1 '} as input excitations at 
flight condition i. 

Define vectors {y/, y/\ . . . , y/ 7i } as response measure - 
10 ments at flight condition i. 

Define matrices {P,, P 2 , . . . , P„,} as plant models at flight 
condition i. 

Choose initial A, 


for j = 1 : n f { 

for i = 1 : n* { 

Validate F u (P,. A) using y\ and u\ 
increase size of A if necessary to validate 

} 

} 

for i=l : n f { 

compute flutter margin u(F u (P,, A)) 


flight condition i. 

Define vectors {y/, y, 2 y/ 1 '} as response measure- 

ments at flight condition i. 

Define matrices {P,, P 2 as plant models at flight 

condition i. 


for i = l : rif { 

Choose initial Aj 
for j = 1 : n s { 

Validate F u (Pj, Ap using y', and \i', 
increase size of A, if necessary to validate 

) . 35 
compute flutter margin S'flutterfrom ti(F u (P h Ap) 

} 

A global approach uses the entire set of flight data from 
test points throughout the flight envelope to generate a single 40 
uncertainty description for all nominal aircraft models. The 
magnitudes of the uncertainty operators are chosen such that 
all nominal models with the associated uncertainty descrip- 
tion are not invalidated by any of the flight data sets. 

Several advantages and disadvantages exist to using the 45 
global approach instead of the local approach. One disad- 
vantage is a possible large increase in conservatism of the 
flutter margin because the uncertainty description is not 
minimized at each flight condition. A single particularly 
inaccurately plant model will require large uncertainty 50 
operators that may be highly conservative for plant models 
at flight conditions that are better representations of the true 
dynamics. 

One advantage to this approach, however, is that the 
uncertainty description is truly worst-case with respect to the 55 
entire flight envelope. The worst-case errors from the worst- 
case fight condition are used to generate the uncertainty 
description for all conditions. Also, this approach is not very 
sensitive to poorly measured flight data. A poorly modeled 
modal response may only appear in certain data sets. The 60 
local approach would not include uncertainty for these 
dynamics at conditions that did not clearly observe this 
modal response, so the resulting flutter margin would not 
account for the true level of modeling errors. The flutter 
margin generated with the global approach may be more 65 
conservative than the local approach, but this approach 
introduces a corresponding higher margin of safety. 
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Hybrid approaches have also been formulated that mix the 
local and global approaches. One straightforward hybrid 
approach is to generate an uncertainty description using all 
data from a small range of flight conditions. This approach 
can be useful for separately considering sets of plant models 
that are generated using different techniques. For example, 
the model-generating package used for this paper computes 
all subsonic plant models with a doublet-lattice algorithm 
and the supersonic models are generated with constant-panel 
algorithms. A hybrid approach could be used to reflect this 
knowledge and consider groups of subsonic, supersonic, and 
transonic plants independently. 

The approaches outlined here are certainly not exhaustive. 
A weighted-norm approach can be formulated that uses 
flight data from the entire envelope but depends heavily on 
a particular subset of that data. Other approaches could 
concentrate on particular dynamics through modal filtering 
techniques to generate separate uncertainty descriptions for 
individual modes. 

Although particular embodiments of the invention have 
been described and illustrated herein, it is recognized that 
modifications may readily occur to those skilled in the art. 
For example, although matrix algebra has been used 
throughout to describe the invention in order to facilitate its 
implementation by computer programming since the matri- 
ces define for the programmer how the elements are to be 
processed, the invention is not limited to the implementation 
disclosed in this manner. Any other form of mathematical 
expression and its implementation by computer program- 
ming will be equivalent. Consequently, it is intended that the 
claims be interpreted to cover such modifications and 
equivalents thereof. 

What is claimed is: 

1. An on-line method for robust flutter prediction in 
expanding a safe flight envelope for an aircraft to define an 
envelope of safe flight conditions, comprising two parts, a 
first part comprising the steps of 

(al) generating a computer model of said aircraft, 

(a2) computing flutter margins of said aircraft from said 
computer model using a singular value \i to define a first 
safe flight envelope and optionally using a well known 
traditional method for defining a second safe flight 
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envelope used in double checking the safety of said first 
safe flight envelope defined by its computed flutter 
margins, 

(a3) if computed flutter margins thus double checked are 
found to define an envelope of safe flight conditions, or 
an optional method of defining a safe flight envelope is 
not used, proceed to a second of said two parts, 
said second part comprising in-flight steps of, 

(b 1 ) taking said aircraft to a safest point within said first 
safe flight envelope at a flight condition, F, and mea- 
suring flight data at said present condition F, including 
dynamic pressure defined by parameters comprising 
altitude and air speed; 

(b2) compute a predicted flutter point, F ; ,, at a higher 
dynamic pressure than at said condition F using an 
algorithm based on said flight data and said singular 
value p; 

(b3) determine dynamic pressure difference between 
dynamic pressure of said aircraft at said present flight 
condition F and at said predicted flutter point F /f \ 
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(b4) if said dynamic pressure difference is large, take said 
aircraft from present flight test condition to a new flight 
test condition F^F+A^, where A F are the changes to 
parameters of flight condition F comprising altitude to 
take said aircraft to a new flight condition F p from flight 
condition F, to incrementally expand testing said first 
safe flight envelope and repeat steps b 1 , b2, b3 and b4; 
but if said dynamic pressure difference is small, declare 
to the last predicted flutter margin F^ as a point on said 
aircraft’s expanded safe flight envelope, and repeat 
steps bl-b4 until sufficient flight conditions within said 
first safe flight envelope have been tested to expand 
margins thereof; 

15 

whereby a robust expanded flight envelope is determined for 
said aircraft using said singular value p and said in-flight 
steps. 

* * * * * 



